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ABSTRACT 


The Marine Corps uses its Selective Reenlisted Bonus (SRB) Program to 
influence Marines to reenlist for a designated term into certain Military Occupational 
Specialties (MOSs) in order to reach planned manpower goals. The bonus amount is 
determined by selecting an “SRB multiplier” for each combination of MOS and Zone 
(“MOSZ”). (“Zone” corresponds to length of service.) A higher multiplier means a 
larger bonus and leads to a higher percentage of Marines reenlisting. That percentage, 
predicted by an existing forecasting model, is assumed exact here. 

The “SRB multiplier model” assigns multipliers to minimize a sum of weighted 
squared deviations from MOSZ targets subject to a budget constraint. This model is 
implemented as a generalized assignment problem, and solved approximately on a 
personal computer using Lagrangian relaxation and a secondary heuristic. (The 
algorithm is programmed in Visual Basic for Applications and has an Excel interface.) 

Data for FY04 shows 491 bonus-eligible MOSZs. With up to 11 possible 
multiplier values, this yields a model with 5,401 0-1 variables and 491 constraints. A 
solution within 0.0018% of optimality is reached in 1.4 seconds on 1.58 GHz personal 
computer. Standard integer-programming software verifies the correctness of the 
solution. 
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EXECUTIVE SUMMARY 


The United States Marine Corps (USMC) uses its Selective Reenlisted Bonus 
(SRB) Program to influence Marines to reenlist for a designated term into certain 
Military Occupational Specialties (MOS) in order to reach planned manpower goals. The 
amount of the bonus is determined by selecting an SRB multiplier for each combination 
of MOS and Zone (“MOSZ”). (“Zone” partitions Marines into three groups based on 
current length of service.) A higher multiplier means a larger bonus and leads to a higher 
percentage of eligible Marines reenlisting. That percentage is predicted by an existing 
forecasting model, whose predictions are assumed exact. 

The total amount of SRB funding available for a fiscal year (FY) is limited by the 
amount approved and allocated in the FY Defense Budget. Thus, not all MOSZs can be 
assigned the multiplier necessary to reach reenlistment goals. Hence, it is important to 
allocate the available SRB funding to those MOSZs that are critical to accomplishing the 
Marine Corps’ mission. So, the problem that is addressed by this thesis is to determine 
the optimal set of multipliers to reach reenlistment goals. 

The “SRB multiplier model” assigns multipliers to minimize a sum of weighted 
squared deviations from MOSZ targets subject to a budget constraint. This was described 
and implemented as a generalized assignment problem in a 1986 Naval Postgraduate 
School thesis by U.S. Marine Corps Captain Dean D. DeWolfe, entitled Determination of 
Selective Reenlistment Bonus Multipliers in the United States Marine Corps. This model 
has not been used by the USMC, however, because it is no longer current with respect to 
SRB rules and parameter definitions, and because it was not implemented in a user- 
friendly programming environment. 

This thesis overcomes the earlier problems with the SRB multiplier model by (a) 
updating parameter definitions and other model constructs to represent the current SRB 
environment, and (b) creating a flexible implementation of DeWolfe’s Lagrangian-based 
solution algorithm on a personal computer. The flexible implementation should allow the 
model and solution algorithm to track changes in data and structure for the SRB program 
as it evolves in the future. The algorithm is programmed in Visual Basic for Applications 
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and is supplied with a user interface through Microsoft Excel. The algorithm first 
optimizes a Lagrangian function that relaxes the budget constraint, and saves the best 
feasible solution it finds during this process. The algorithm then attempts to improve that 
solution with a greedy marginal-rate-of-return heuristic. 

Data for FY04 yields 491 MOSZs that are eligible for an SRB multiplier. With 
up to 11 multiplier values for some MOSZs, this results in a model with 5,401 0-1 
variables and 491 constraints. A solution is reached in 1.4 seconds on a 1.58 GHz 
personal computer. The solution is within 0.0018% of optimality and consumes 99.99% 
of the budget. 
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I. 


INTRODUCTION 


Each year the United States Marine Corps (USMC) must set the Selective 
Reenlistment Bonus (SRB) for each Military Occupational Specialty (MOS) and Service 
Time Zone (Zone) (USMC 1990). The goal is to selectively encourage Marines to 
reenlist and enable the USMC to reach reenlistment targets for the given fiscal year. The 
purpose of this thesis is to give manpower planners at Manpower & Reserve Affairs 
(M&RA), Headquarters Marine Corps, a tool that will identify the optimal “SRB 
multiplier,” for each MOS and Zone combination (MOSZ). Theses multipliers set the 
monetary values of the bonuses. 

A. ENLISTED MARINES, MOS AND ZONE 

The Marine Corps uses two models, the First Term Alignment Plan (FTAP) and 
the Subsequent Term Alignment Plan (STAP), to determine the number of reenlistments 
that are necessary to meet the future force structure as depicted by the Grade Adjusted 
Recapitulation (GAR) (Manpower Plans and Policy Branch-Integration and Analysis 
2007). The GAR is a planning tool that specifies target enlistment numbers for all MOSs 
and ranks. For planning purposes, a Marine is classified by pay grade and time in 
service. Pay grade is associated with a Marine’s level of responsibility while time in 
service correlates to experience. A Private (E-l) is the lowest enlisted pay grade, is the 
least experienced Marine with less than one year in service, and is generally still in initial 
training. A Sergeant Major (E-9) is the highest enlisted pay grade and is the most 
experienced enlisted Marine, having 15 or more years in service. 

For reenlistments purposes, manpower planners are concerned with Marines who 
are in or about to enter the career force, meaning they have 21 or more months time in 
service. These Marines are placed into one of three Service Time Zones according to 
time in service as listed in Table 1. The FTAP model covers Marines in Zone A while 
the STAP model covers Marines in Zones B and C. 
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Table 1. Service Time Zones designated by time in service. 


ZONE 

FROM 

TO 

A 

21 Months 

6 Years 

B 

6 Years 

10 Years 

C 

10 Years 

14 Years 


FTAP is the retention model used to help determine the number of reenlistments 
needed from Zone-A Marines, according to MOS, to meet the career-force requirements 
set by the GAR (USMC 2004). A Zone-A Marine is in his initial enlistment contract, 
including extensions. The purpose of the FTAP is to fill vacancies in the career force that 
result from separation or loss for any reason while preventing promotion stagnation and 
ensuring opportunities for advancement. 

STAP is similar to the FTAP, but focuses on subsequent reenlistments of Marines 
who constitute the career force (Manpower Plans and Policy Branch-Integration and 
Analysis, 2007). The primary purpose of the STAP is to “move” career-force inventory 
levels towards required levels. STAP recognizes that any career-force imbalances cannot 
be “fixed” in one year, and tries to make the following year better in relation to the 
imbalance in the previous year. 

FTAP is important because it covers the most flexible opportunity the Marine 
Corps has to shape the career force, with respect to the number of Marines in each MOS 
and pay grade. Only during a Marine’s first reenlistment opportunity can the Marine 
Corps deny reenlistment (without monetary cost) because his particular MOS is full: A 
reenlistment is allowed only when there is a specific vacancy in the career force that 
needs to be filled. 

Provided that a first-term Marine is qualified, any Marine who is denied 
reenlistment in his original MOS, because no vacancies exist, is allowed to reenlist in 
another MOS that needs additional manpower (USMC 1993). This process is called a 
Lateral Move or “Lat Move.” Reenlistments for subsequent-term (Zone-B and Zone-C) 
Marines are not based on vacancies but instead are based on targets. Furthermore, the 
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MOS will not close to reenlistments even if the target is reached. Although allowed in 
specific cases, Lat Moves are not used with subsequent-term Marines as they are with 
first-term Marines. 

The ability to deny first-term reenlistment in conjunction with the use of Lat 
Moves allows manpower planners to control the number of Marines that advance into the 
career force from the first-term population, in specific MOSs, and to maintain the 
personnel goals set forth by the GAR. During subsequent reenlistment opportunities with 
a career-force Marine, the Marine cannot be denied reenlistment provided the Marine has 
no performance or criminal-conduct problems and has been promoted on schedule. If, in 
order to control the MOS population, manpower planners need to separate a subsequent- 
term Marine who has no difficulties, the affected Marine is entitled to a substantial 
involuntary- separation allowance. 

Historically, most MOSZs do not need a bonus to encourage enough Marines to 
reenlist to meet assigned targets. Approximately half of the bonus eligible MOSZs will 
not receive a bonus because the anticipated manpower shortfall is not severe enough to 
warrant a bonus. For the remaining MOSZs, the manpower shortfall is critical and the 
only option that can help meet the remaining reenlistment requirements is the SRB. 

B. HOW THE SRB IS USED 

When too few Marines want to reenlist or Lat Move into a critical MOSZ, the 
Marine Corps uses the SRB to increase the number of reenlistments there. Typically, a 
critical MOS falls into one of the following categories (Morgan 2006): 

• Technical skills requiring large training and/or replacement costs, 

• Skills in high demand in the civilian sector, 

• Skills for which recruiting is difficult, 

• Skills that are crucial to combat readiness, and 

• Skills in high demand but with a small reenlistment population. 
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The Marine Corps separates Marines into one of three Service Time Zones by 
their time in service, as defined in the previous section. An SRB can be awarded to a 
Marine only once in each Zone, but a Marine may be eligible to receive an SRB in each 
Zone through his career. 

The monetary amount awarded for an SRB is computed as a function of the 
Marine’s pay, additional years of service that reenlistment obligates, and the “strength” of 
the bonus referred to as an SRB multiplier. More precisely, the SRB is defined as 

The Minimum of: 

Monthly Base Pay x Years of Additional Obligated Service 
x SRB Multiplier for MOSZ 

and 

Maximum Allowable Bonus for the Fiscal Year. 


The Marine Corps uses SRB multipliers ranging from 0 to 5 in increments of 0.5. 
In an MOS that does not need an incentive; the multiplier is 0 and no bonus is offered. In 
an MOS with an expected shortfall, the multiplier will typically be non-zero. 

For Zone A, manpower planners have the power to control the amount of total 
bonus money that is paid by closing an MOS to reenlistment. Manpower planners cannot, 
however, close off bonuses or reenlistments to a Zone-B or C MOS, once a bonus is 
offered there. This means that if all of the eligible Marines in an MOS reenlist, all will 
receive the bonus money. These reasons combined with the limited yearly budget for the 
SRB Program make picking the appropriate multiplier values crucial. 


C. CURRENT SRB MODELS 


Currently, the SRB Program uses two manual heuristics to specify SRB multiplier 
values, one heuristic for FTAP (Zone A) and the other for STAP (Zones B and C). The 
total budget is first divided between FTAP and STAP, and then the heuristics are run 
separately. Both heuristics are implemented using Excel spreadsheets that first set the 
multiplier in each MOSZ at the minimum level that is estimated to lead to the meeting of 
its reenlistment requirement. In order to determine which multiplier to assign, both 
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models use expected reenlistment rates provided by the Center for Naval Analysis (CNA) 
(North 1995). Once the initial settings are complete, the expected costs for the selected 
multipliers are summed over all MOSZs to estimate a total cost of the program. If this 
estimate exceeds the SRB budget, individual SRB multipliers for less critical MOSs are 
manually reduced until it does not. 

Manpower planners combine the output of FTAP and STAP into a list of 
prospective SRB multipliers, and carry out an iterative, negotiating process with other 
departments at Headquarters Marine Corps, to agree upon the final SRB multiplier 
values. The other departments include the occupational-field sponsors, assignment 
monitors, command-level representatives, and career-retention specialists, all of whom 
are the end-users of the SRB Program. 

D. SUMMARY OF PAST WORK ON SRB OPTIMIZATION AND OTHER 

RELATED RESEARCH 

A 1986 Naval Postgraduate School thesis by U.S. Marine Corps Captain Dean D. 
DeWolfe, entitled Determination of Selective Reenlistment Bonus Multipliers in the 
United States Marine Corps , is the basis for this thesis. First, DeWolfe sets up and 
defines the problem of selecting an SRB multiplier for each MOSZ as a nonlinear 
knapsack problem. The objective is to select multipliers that minimize a total weighted 
deviation from the number of estimated reenlistments versus the reenlistment targets, 
subject to a budget constraint; all coefficients are expressed as expected values. Next, 
DeWolfe reformulates the problem as a special generalized assignment problem shows 
how to maximize a Lagrangian lower bound on the optimal objective value and, in the 
process, identify a good feasible solution. Finally, due to the nature of the Lagrangian 
relaxation, it is possible the solution does not use the entire budget. In that case a 
secondary heuristic tries to increase multipliers to consume it all. The heuristic uses a 
marginal rate of return to determine which MOS multiplier value can be raised to the next 
level to achieve the largest decrease in the objective value per unit of (additional) budget 
consumed. 
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Although M&RA liked how the model solved the SRB multiplier problem, there 
were problems that made it unusable. First, the model was written in FORTRAN 77 
which is not familiar to personnel at M&RA. Second, there have been many changes to 
the SRB program rules since 1986 that needed to be updated (SRB multiplier increments, 
budgets, lump sums, maximum allowed bonuses), and many of these rules were “hard¬ 
wired” into the code. Finally, and most importantly, M&RA wanted a model that could 
mn in a familiar computational environment like Excel’s VBA (Visual Basic for 
Applications). Such an environment would allow them to internally use and maintain the 
model. 

Other related research includes two NPS theses by U.S. Marine Corps Captain 
Jonathan D. Raymond (2006) entitled Determining the Number of Reenlistments 
Necessary to Satisfy Future Force Requirements, and by U.S. Marine Corps Major Dean 
G. Conatser (2006) entitled Forecasting U.S. Marine Corps Reenlistments by Military 
Occupational Specialty and Grade. Raymond created a model to determine the number 
of reenlistments that are necessary for each MOS and grade by applying transition rates 
to the current inventory of Marines who are not eligible for reenlistment and subtracting 
that amount from the desired force structure in the GAR. Conatser created a model to 
predict the number of reenlistments for both FTAP and STAP Marines by MOS and 
grade. These two models, developed at the request of M&RA, are the beginning of the 
development of the Career Force Retention Model (CFRM) which is being explored to 
replace the FTAP and STAP models. The CFRM, in conjunction with this thesis, will 
help manpower planners estimate the required number of new recruits, reenlistments, and 
the most effective way to set SRB multipliers to influence the reenlistments. 

This thesis extends and updates DeWolfe’s model and solution algorithm for use 
by M&RA. This is accomplished by revising budget constraints, adding a Lat Move 
functionality, and making all model parameters accessible to the user (i.e., no hard-wiring 
of parameter values into the code). The ultimate goal of this thesis is to create a model 
that is user-friendly and “user-updatable.” 

In a paper by DeWolfe, Stevens, and Wood in 1993 entitled Setting Military 

Reenlistment Bonus, the authors expand upon the model that was developed in 
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DeWolfe’s 1986 thesis for use by the United States Army. First a more complicated 
formulation was created by adding an additional constraint that limited the number of 
“high-level” bonuses to 10% of the total bonuses offered. Second, the formulation found 
multipliers for each MOS and Zone based on expected reenlistment rates for 3, 4, 5, and 
6 years of additional obligated service. Since the 10% constraint is no longer current and 
the Marine Corps only offers four-year reenlistments, the more complicated formulation 
is not necessary in this thesis. 

E. THESIS OUTLINE 

In Chapter II, the SRB multiplier problem is formulated as a generalized 
assignment problem. The data and variables are described along with the methodology 
used to compute model coefficients like the training costs in the object function. This 
model is then implemented and solved in the General Algebraic Modeling System 
(“GAMS”; see Rosenthal 2007), with an Excel interface, for testing purposes. Although 
the GAMS implementation and solution is a perfectly reasonable one, it is unsatisfactory 
for the final product of this thesis, because it relies on an unfamiliar computational 
environment for manpower planners. Instead, this solution is used in Chapter III to check 
the accuracy of a solution method based on Lagrangian relaxation, whose implementation 
and testing in VBA and Excel is also described in Chapter III. Chapter IV presents 
conclusions and recommendations for further work. Appendices A and B contain the 
source code for the GAMS and Excel VBA models, respectively. Appendix C contains 
sample inputs and outputs from the Excel VBA model. 
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II. GENERALIZED ASSIGNMENT MODEL 


This chapter models the optimization of SRB multipliers across all MOSZs as a 
generalized assignment problem (GAP). First, the model terms and coefficients are 
defined. Second, the differences from DeWolfe’s model and the current model are 
explained. Finally, the estimation methods used for some of the coefficients in the model 
are defined. 

A. PROBLEM FORMULATION 

There is a total of m MOSZs that encompass all MOSs and Zones. Each MOSZ 
will take on one of n multiplier “levels,” with level j = 1 corresponding to a multiplier 
value of Vj = 0. Current SRB policy uses multipliers in increments of 0.5, up to a 
maximum of 5.0, for a total of 11 multiplier levels. Thus, each increase in level j results 
in an increase of 0.5 in v ; . For example, level j = 2 corresponds to v 2 = 0.5 , j = 3 to 

v 3 =1.0, and the maximum level of j = 11 corresponds to v u = 5. The model is not 
limited to a predefined number of multiplier levels, as long as v . +1 - v . is constant for all 
levels j. 

Currently, there is no set maximum multiplier value assigned to any single MOSZ 
or group of MOSZs that is different from the overall maximum multiplier mentioned 
above. However, the multiplier cannot exceed a value that would result in a bonus that 
exceeds the maximum allowable bonus. The maximum allowable multiplier is easy to 
compute and is represented through n i , the maximum multiplier level for MOSZ i. 

Each MOSZ has an associated cost (budget consumption) and a penalty 
(objective-function value) that depends on the expected number of reenlistments for the 
chosen multiplier. The overall objective is to assign each MOSZ a multiplier that 
minimizes the total penalty and keeps the total cost under budget. The following 
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formulation of the SRB multiplier problem as a GAP is similar to DeWolfe’s but some 
variables and parameters have been changed to account for Lat Moves and the use of 
non-integer multipliers. 

1. Indices 

i = MOSZ 

j = 1allowable SRB multiplier levels for MOS i 

2 . Parameters [units] 

d tj reenlistment-target deviation penalty for MOSZ i with 

multiplier level j [arbitrary penalty units] 

budget outlay incurred by setting multiplier level to j for MOSZ i [dollars] 

B budget [dollars] 

3. Decision Variables 

Xy 1 if multiplier level j is selected for MOSZ i, 0 otherwise 

4. Formulation for GAP 

m >’/ 

z*=min Y.Y.dijXij ( pl ) 

i= 1 7=1 

m n , 

'• L Eiw 58 («) 

1=1 7=1 

E^,=l Vi (C2) 

7=1 

x ij e {0,1} Vz, j (C3) 

The objective function (PI) sums penalties across all MOSZs. The penalty is 0 
for a reenlistment target that is met exactly, so the best possible objective-function value 
is 0. The first constraint (Cl) limits total bonus expenditures to the available budget (at 
least in terms of expected values). Note that c n =0 for all MOSZs i, that is, no cost is 
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incurred if no bonus is offered. The second constraint (C2) ensures that each MOSZ 
receives exactly one multiplier. This model follows DeWolfe’s exactly, except that 
coefficients d and c are computed differently, as described below: 


W, r 

d = MTC X 


max {0,(D. - LM j LMF j ) - RRR ij E i ! 

+ Qx max {o, RRR, E i - (Z> - LM i LMF i )} 2 


(Eqn 2.1) 


c tj = RRR tj E^VjY 


(Eqn 2.2) 


where coefficients are defined by 


D j number of Marines needed to reenlist in MOSZ i [Marines] 

IT exogenous weighting factor for MOSZ i [unitless] 

I] average monthly pay for Marines in MOSZ i [dollars] 
vj numerical multiplier value for multiplier level j 

RRR j expected reenlistment fraction for Marines in MOSZ i and multiplier level j [fraction] 
B budget [dollars] 

Q relative penalty weight of overages to underages [fraction] 

Y average years of reenlistment [years] 

MTC maximum training cost for all MOSZs [dollars] 

LM j number of MOSZ i school seats open for Lateral Moves [school seats] 

LMF i fraction of available school seats to use for Lateral Moves into MOSZ i [fraction] 

7] training cost of MOSZ i [dollars] 

Aj total number of Marines in MOSZ i [Marines] 

E i number of Marines eligible for reenlistment in MOSZ i [Marines] 


In Eqn 2.1, a Lat-Move functionality is added by the term LM j LMF j . This term 
can reduce the number of reenlistments that is needed from the eligible group of Marines 
in each MOSZ. Each MOSZ has a predetermined number of MOS school seats, LM t , 

that are open for Lat Moves. Depending on the MOSZ, the number of school seats can 
range from zero to the total number of Marines that are needed to reenlist. In some 
MOSZs, manpower planners will want to use all available school seats and set 
LMF i =1.0. In others, planners will want more control over the number of new Marines 
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that LatMove into the MOSZ and set LMF j < 1. Manpower planners can control the 
number of LatMoves by assigning LMF n a fraction of the school seats to use for each 

MOSZ. Historically there are always enough Lat-Move eligible Marines to fill school 
seats in a given FY, but manpower planners will need to verify that this assumption is 
still valid. 

There are two changes to Eqn 2.2 over DeWolfe’s formulation. First, the total 
bonus amount for each MOSZ is used in the budget constraints to account for the use of 
lump-sum payments instead of the anniversary payments in use in 1986. Second, a new 
term v. is introduced to account for the use of non-integer SRB multipliers. It is a 

parameter for the SRB multiplier value that is used in calculating the expected cost of 
using an arbitrary multiplier v.. With this new parameter, manpower planners will not 

be limited to the integer-only multipliers of 1986 or the multipliers increments of 0.5 in 
use today. 

B. TRAINING-COST DATA 

A key term in the weighting function for the objective of the SRB model is the 
training-cost parameter, which weights the deviation penalty linearly with increasing 
training cost: If a shortfall of X Marines must occur in some MOSZ, it is better for that 
to occur in an MOS where the “replacement cost,” based on training costs, is lower. 
Unfortunately, the Marine Corps no longer tracks training-school costs, because multiple 
sources of funding makes this accounting task difficult. Currently, the Marine Corps can 
only supply the school costs for non-Marine Corps schools because the Marine Corps 
pays a known amount for each Marine to attend. In order to overcome the lack of actual 
cost data, a new training-cost estimate is developed here, based on Zone, length of school 
or schools attended, graduation rate, and an estimate of the student’s daily pay. All 
MOSZ training-cost parameter estimates use a variation of the following equation: 

Training-cost estimate = 

(Training Days x Daily Pay) / Graduation Rate (Eqn 2.3) 
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1 . 


Training-Cost Estimate for Zone A 


Zone-A training costs are broken into two categories: Initial training MOSs and 
Lat-Move-only MOSs. Both categories use Eqn 2.3 to estimate corresponding training 
costs but use different daily pay amounts to reflect the average pay received by a Marine 
while in training. This thesis applies, for initial training MOSs, the daily pay for an E-l 
with less than two years in service; for Lat-Move-only MOSs, it uses the daily pay for an 
E-5 with more than four years in service. 

2. Training-Cost Estimate for Zones B and C 

Zone-B and -C training costs are also based on Eqn 2.3 but are more complicated 
because of the various tracks that a Marines’ career can take in his MOS, an MOS that 
may, in fact, change. First, the Marine can be promoted and keep his original MOS. 
Under those circumstances, training cost is based on the MOS’s Zone-A training cost. 
Second, a Marine can be promoted and receive a new MOS to distinguish new 
responsibilities. The new MOS can be supplied by a single MOS or by multiple “feeder 
MOSs.” In the former case, the training-cost estimate is also based on the Zone-A 
training cost of the old MOS; otherwise the training cost is the average total training costs 
of the Zone-A feeder MOSs. In all of the above cases, it is possible that the Marine will 
be required to attend another MOS training school. If applicable, the new training cost 
for the school is estimated through Eqn 2.3 and is added to the prior-training cost 
estimate that is determined by the MOS track. 

C. GAMS MODEL 

Using FY2004 data, the GAP is implemented and solved in GAMS using the XA 
solver (GAMS 2007). The FY2004 data set contains all 809 USMC MOSZs, of which 
491 MOSZs need a multiplier greater than zero to reach their reenlistment goals. The 
final objective value is 57.8525 is obtained in 0.02 seconds of solver time, and is proven 
to be within 0.0006% of optimality. Of the 491 MOSZs, 385 receive a multiplier greater 
than zero. Table 2 shows the distribution of the eligible MOSZs and their selected 
multipliers. 
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Table 2. Distribution of SRB Multipliers in GAMS Solution 


Multiplier 

O 

o 

0.5 

1.0 

1.5 

2.0 

2.5 

o 

CO 

LO 

CO 

4.0 

4.5 

o 

LO 

MOSZ 

106 

39 

99 

41 

33 

95 

74 

0 

1 

3 

0 


The eligible MOSZs have a reenlistment goal of 8,887 Marines, and with the 
chosen multipliers 5,620.75 Marines are expected to reenlist. The SRB budget for 
FY2004 is $54 million and only $474 goes unspent in the GAMS solution. 
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III. SOLVING THE GENERALIZED ASSIGNMENT MODEL 
WITH LAGRANGIAN RELAXATION 


Although the updated GAP formulation of the SRB multiplier problem solves 
quickly in GAMS, manpower planners cannot readily use this implementation due to lack 
of access to GAMS. In order to supply manpower planners with a model they can use, 
the GAP model is implemented and solved here using the Lagrangian-relaxation 
technique described by DeWolfe, but using a current programming language that can mn 
in a computational environment that manpower planners have access to. The following 
section borrows heavily from DeWolfe (1986), and the reader should see that thesis for 
more details. 

A. LAGRANGIAN RELAXATION 

The budget constraint (Cl) in the GAP can be viewed as a “complicating 
constraint.” That is, the model solves easily if (Cl) is relaxed. Using Lagrangian 
relaxation, the budget constraint can be moved to the objective function using a 
Lagrangian multiplier that acts as a penalty for constraint violation, or as an incentive to 
spend. The relaxed formulation is: 

For A > 0, 

m a j m a j 

LR(2): z(A) = minY J Y J d ij x ij + ACZ^c^-B) (P3) 

M j =1 i'=l j =1 

n i 

s.t. v * 

x u e { 0 ,!} 7 

The inner portion of P3 can be simplified and rewritten as: 

m rij 

min^] J ( d u + ^ c ij ) x ij - ^ B ( E 9 n 3.1) 

M j =1 
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Thus, for fixed X , z(X) is computed by solving a simple selection problem: For 
each MOSZ i , choose x. to be 1 for the largest value of cL + Xc tj , and set other x. to 0. 

It is well known, then, that z(X) < z* , i.e., z(X) provides a lower bound on z*. 

The best lower bound is found by maximizing z(X) with respect to X to obtain 
X* . This can be done by bracketing the optimal X in some range [L t , L u ], and then 
narrowing the bracket through bisection search. L t = 0 is always valid. For X 
sufficiently large, the solution to LR(/l) will spend as few SRB budget dollars as 
possible. This will occur when setting x, j = 1 for all i is the optimal action, and this will 
occur when 

d n + Xc n < d tj + Xc tj Vi, / > 2 (Eqn 3.2) 

Since c n = 0 , this leads to 

d n <d ij +Xc ij Vi,j>2 

Thus 

4 = max {(d n - d tj ) / Cy } V/, j > 2 

is a valid value for the right-hand side of the bracket. 

The bisection search is based on the fact that if the budget constraint is violated 
for a given X , then X* > X ; otherwise, X* < X . 

The above methodology is implemented just as in DeWolfe (1986), but using 
VBA, and the best Lagrangian lower bound, z(X*) , is established. An upper bound on 
the optimal solution to the SRB multiplier problem is found by examining the feasible 
solution sets to z(X) that are encountered in the course of maximizing z(X). (When 

X > X* , the solution to LR(Z) is guaranteed to be feasible with respect to the budget 
constraint.) 


(Eqn 3.3) 

(Eqn 3.4) 
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The best solution identified while optimizing z(Z) may not consume the entire 
budget. If unspent budget remains, it may be possible to improve the solution using the 
marginal-rate-of-retum heuristic described by DeWolfe. The marginal rate of return for 
MOSZ i is defined as: 

ROR t = (d„ -d, J+1 )/(c iJ+l -Cy) (Eqn 3.5) 

where the current feasible solution hasxy =1. The numerator in Eqn 3.5 calculates the 
amount z(Z) can be improved by increasing the SRB level for MOSZ i from its current 
level of j to j + 1; the denominator calculates the corresponding additional budget 
expenditure. In essence, ROR j calculates the “bang for the buck” for raising the 

multiplier level for each MOSZ individually. The heuristic repeatedly applies these rules 
until no additional budget can be feasibly consumed: 

a. Given the current values for x i} , compute ROR i for each MOSZ i. 

b. Select MOSZ i* with the largest value of ROR i such that (a) setting 
Xy = 0 and x t /+1 = 1 does not violate the budget constraint, and (b) j +1 
does not exceed the largest allowable SRB level for the MOSZ. 

c. Set Xy =0 and x t j+i =1. 

B. EXCEL VBA MODEL SOLUTION BY LAGRANGIAN RELAXATION 

Using FY2004 data, the Lagrangian-based solution method for the GAP, as 
described above, is implemented in EXCEL using VBA. For the FY2004 data described 
previously, the final objective value of 57.8691 is obtained in 1.4 seconds, and is proven 
to be within 0.0018% of optimality. Of the 491 MOSZs eligible for a bonus, 385 receive 
a multiplier greater than zero. Table 2 shows the distribution of the eligible MOSZs and 
their selected multipliers. 

Table 3. Distribution of SRB Multipliers in VBA Solution 


Multiplier 

O 

o 

LO 

o 

1.0 

1.5 

2.0 

2.5 

o 

CO 

3.5 

4.0 

4.5 

o 

LO 

MOSZ 

106 

39 

94 

42 

35 

97 

74 

0 

1 

3 

0 
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The eligible MOSZs have a reenlistment goal of 8,887 Marines, and 5,621.51 
Marines are expected to reenlist with the chosen multipliers. Of the $54 million budget, 
only $574 goes unspent. 

C. COMPARISON OF DIRECT AND LAGRANGIAN SOLUTIONS 

The VBA and GAMS solution methods produce very similar outputs, with the 
VBA model having the higher objective-function value by 0.0166 units, or by 0.0287%. 
Of the 491 MOSZs that are eligible for an SRB, both models assign non-zero multipliers 
to the same 385 MOSZs. Of these MOSZs, the VBA solution assigns only 11 multiplier 
values that differ from the GAMS solution. Of these, nine multipliers increase by 0.5 in 
the VBA solution, one multiplier decreases by 0.5, and one multiplier increases by 1.5. 

The largest difference in objective-function values results from the lone 
multiplier that decreases by 0.5. The decrease raises the objective value by 0.1057 units 
or 0.18% of the total VBA objective-function value. For the MOSZ with the largest 
multiplier difference in the VBA solution, the multiplier is raised from the GAMS- 
solution value of 1.5 to 3.0. Despite this large difference, the impact on the objective 
value is negligible: The objective value decreases by 0.1106 units or 0.06% of the total 
objective value. This small change in the objective value results from the fact that the 
relevant MOSZ contains only two Marines, and so even the worst multiplier assignment 
contributes little to the overall objective-function value (penalty). 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


In summary, the model developed in this thesis provides Marine Corps manpower 
planners a new tool to help manage the Selective Reenlistment Bonus (SRB) program 
(USMC 1990). The model’s solution specifies “SRB multipliers” for each Military 
Occupational Specialty (MOS) and Service Time Zone (Zone) combination that define 
the bonus amount, and are projected to encourage reenlistments that will (a) meet annual 
reenlistment targets as best possible, and (b) not exceed the available SRB budget. The 
new extended SRB optimization model and solution procedure incorporates key aspects 
of the methodology developed in DeWolfe (1986), but makes these improvements: 

• Gives manpower planners the ability to exploit Lateral Moves, 

• Allows a more flexible set of multiplier values that need not be integers, 

• Makes budget and other parameters “user updateable,” and 

• Implements the model and solution procedure in a familiar computing 
environment for manpower planners, and at no cost (beyond software and 
hardware that these planners already possess). 

The new model does have its limitations. When using the Lateral-Move 
functionality, the critical assumption about the availability of Marines eligible for “Lat 
Moves” needs to be checked. Also, the expected number of reenlistments for a given 
multiplier is based on a forecasting model that is assumed perfect. 

Two areas that could improve the effectiveness of the model warrant further 

study: 

• The development of a database that more accurately reflects true training 
costs, rather than the one used in this thesis, which estimates those costs 
based on training days and average pay, 
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• The expansion of the model to better forecast the expected cost of a given 
multiplier by breaking down the Marines in a MOSZ according to pay 
grade. 
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APPENDIX A: GAMS SOURCE CODE 


This appendix contains the source code for the nonlinear optimization problem 
solved in GAMS as a comparison and check for the Excel VBA solution. 

STITLE SRB Multiplier Optimization 

*-DEFAULTS- 

SOFFUPPER OFFSYMLIST OFFSYMXREF 

OPTIONS 
LIMROW = 0 
LIMCOL = 0 
ITERLIM = 500000 
RESLIM = 100000 
SOLPRINT = OFF 
DECIMALS = 2 
LP = XA 
RMIP = XA 
MIP = XA 
OPTCR = 0.001 
OPTCA = 0 


SONTEXT 

Original: 20 Jan 2007 

Author : Kent A. Robbins, Jr. 

Class : Thesis 
Revised: 10 August 2007 

-Deleted Zone Loop (Now just associated with MOS) 

-Added read-in data function 

-Dropped Bonus limit constraint which is accounted for 
in the data that is received from the Excel VBA 
spreadsheet that writes the files. 

-Added fixing variables for Multipliers 
-Added output into Excel file 
-dropped OPTCA from 0.0 to .001 
Description: SRB Multiplier Optimization Program that selects the 
optimum multiplier to minimize the squared weighted 
deviation of the differences between the expected 
and needed number of reenlistments. This program 
is run in shell from Excel's VBA. 
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$OFFTEXT 


*-INDICES - 

SETS 
i MOS 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\MOSZoneList.dat 

m Multiplier 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\MultipleList.dat 


-DATA 


Scalar 


Q "relative weight of overages to shortages wrt reenlistment target" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\Q.dat 
B "Total allowed budget" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\Budget.dat 
Yrs "number of years for reenlistment" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\Yrs.dat 

maxTmCost "maximum training cost of all MOS's" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\MaxTrnCost.dat 


Parameter T(i) "normalized cost of traing a Marine in MOS.Zone(i)" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\TrnCost.dat 


Parameter A(i) "number of Marines in MOS.Zone(i)" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\MOSPop.dat 


Parameter E(i) "number of Marines in MOS.Zone(i) eligible for reenlistment" 
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$INCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\EASPop.dat 


Parameter D(i) "number of Marines required to reenlist in MOS.Zone(i)" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\BoatSpace.dat 


Parameter W(i) "exogenous weighting factor for MOS.Zone(i)" 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\WgtFactor.dat 


Parameter P(i) "average monthly base bay of a Marine in MOS.Zone(i)" 

$INCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\AvgPay.dat 


Parameter MultiplierNumber(m) "corresponding multiplier associated with set (m)" 
SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\MultipleNumbers.dat 


Parameter Preset(i) "gives the value of the preset multiplier for MOS.Zone(i)" 
SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 
Thesis Data Files\PresetMult.dat 


TABLE 

RRR(i,m) "reenlistment response rate for MOS.Zone(i) with Multiplier(m)" 
SONDELIM 

SINCLUDE C:\Documents and Settings\Kent Robbins\My Documents\Thesis\GAMS 

Thesis Data Files\ReenlistmentRates.dat 

SOFFDELIM 

Parameter TEPDev(i,m) "deviation of D(i)-RRR(i,m)*E(i) or 0 (target-expected)"; 
TEPDev(i,m)= max( 0, D(i)-RRR(i,m)*E(i)); 

Parameter ETPDev(i,m) "deviation of RRR(i,m)*E(i)-D(i) or 0 (expected-target)"; 
ETPDev(i,m)= max( 0, RRR(i,m)*E(i)-D(i)); 

Parameter Wght(i) "Weight function for MOS.Zone(i)"; 

Loop((i), IF(A(i) > 0, Wght(i)= (T(i)/maxTrnCost)*W(i)/A(i); 
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ELSE Wght(i)=0;);); 


* -VARIABLE- 

BINARY VARIABLE 

y(i,m) Binary var is 1 if Multiplier m is used for MOS.Zone(i) 
VARIABLE 

Z Total deviation for selection optimal multipliers 

* -EQUATIONS- 

EQUATIONS 

TOTALDEVIATION total cost of SRB allocation 
PickOne(i) only allow one SRB Multiplier 
TotalCost total cost of picking mult m for MOS.Zone(i) 


*-OBJECTIVE FUNCTION- 

TOTALDEVIATION.. 

Z =E= Sum((i,m),(Wght(i)*(TEPDev(i,m)**2+Q*ETPDev(i,m)**2))*y(i,m)) 


*-CONSTRAINTS- 

PickOne(i).. 

Sum(m, y(i,m)) =E= 1; 

TotalCost.. 

Sum((i,m), RRR(i,m)*E(i)*y(i,m)*MultiplierNumber(m)*P(i)*Yrs) =L= B; 


*This loop will fix multipliers if the MOS.Zone(i) as a preset values. 
LOOP((i)$(Preset(i)>-1), 

Loop((m), if(MultiplierNumber(m) = Preset(i), y.fx(i,m) = 1;);); 

); 


MODEL SRB /ALL/; 

SOLVE SRB USING MIP MINIMIZING Z ; 

file out /CADocumcnts and Settings\Kent Robbins\My Documents\Thesis\GAMS Thesis 
Data Files\SRBOptResults.csv/; 
put out; 

*This loop writes the solution to a file that will be read in by Excel. 
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LOOP(i, 

LOOP(m$(y .l(i,m)= 1), 
put MultiplierNumber(m),' 

); 

put /; 

); 


putclose out; 



THIS PAGE INTENTIONALLY LEFT BLANK 
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APPENDIX B: VBA SOURCE CODE 


This appendix contains the source code for the VBA modules that are used to set 

up the SRB optimization problem for both the GAMS and VBA models. This appendix 

also contains the module that solves the Lagrangian relaxation to the Generalized 

Assignment Model solved in GAMS. 

'File: Checks Module in VBA 
'Date Created: 12 December 2006 
'Last Updated: 15 August 2007 

t 

'Checks is a VBA Module that contains a list of subroutines that will perform checks on 
'the data set to verify that it is complete in order to run the optimization software. 

f 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 
Public INumNodes As Long 
Public IMaxIndex As Long 
Public sNodeNames() As String 
Public llndex() As Long 
Public INum E rrors As Long 

'FillZoneABtSpacesO is a subroutine that will compute the correct number of 
'reenlistments for Zone A that are needed when BNA school seats are used. 

f 

'@param - None 

t 

Sub LillZoneABtSpaces() 

Dim rZoneABtSpacesAvail As Range 
Dim rZoneABNASeatsAvail As Range 
Dim rZoneABNASeatUsage As Range 
Dim rZoneABtSpaceAdjust As Range 
Dim iLength As Integer, i As Integer 

Set rZoneABtSpacesAvail = wsMOSZoneA.Range("ZoneABtSpaceStart", _ 
wsMOSZoneA.Range("ZoneABtSpaceStart").End(xlDown)) 

Set rZoneABNASeatsAvail = wsMOSZoneA.Range("ZoneABNASeatStart", _ 
wsMOSZoneA.Range("ZoneABNASeatStart").End(xlDown)) 

Set rZoneABNASeatUsage = wsMOSZoneA.Range("ZoneABNASeatUsageStart", _ 
wsMOSZoneA.Range("ZoneABNASeatUsageStart").End(xlDown)) 
iLength = rZoneABtSpacesAvail.Rows.Count 

Set rZoneABtSpaceAdjust = wsMOSZoneA.Range("ZoneABSAdjStart", _ 

wsMOSZoneA.Range("ZoneABSAdjStart").Offset(iLength, 0)) 
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For i = 1 To iLength 

rZoneABtSpaceAdjust(i) = Application.WorksheetFunction.Max _ 
((rZoneABtSpacesAvail(i) - Application. _ 
WorksheetFunction.RoundUp(rZoneABNASeatsAvail(i) * _ 
rZoneABNASeatUsage(i), 0)), 0) 

Next i 
End Sub 

'MOSZoneCheck() is a subroutine that will verify that the MOS is listed in the Training 
'Cost Worksheet. This subroutine gathers the data and calls subroutine CheckMOS to 
'verify that MOS is listed. If there are MOSs that are not listed a window will pop up 
'and warn the user. 

f 

'@param - None 

f 

Sub MOSZoneCheck() 

'Data Sheet Ranges 

Dim rZoneAMOS As Range 

Dim rZoneBMOS As Range 

Dim rZoneCMOS As Range 

Dim rZoneAMOSTrnVal As Range 

Dim rZoneBMOSTmVal As Range 

Dim rZoneCMOSTmVal As Range 

Training Sheet Ranges 

Dim rTrainMOSAList As Range 

Dim rTrainMOSBCList As Range 

Dim rTrainMOSACostList As Range 

Dim rTrainMOSBCCostList As Range 

'Check & Fill in Zone A Ranges 

Set rZoneAMOS = wsMOSZoneA.Range("ZoneAMOSStart", wsMOSZoneA.Range 
("ZoneAMOSStart").End(xlDown)) 

Set rZoneAMOSTrnVal = wsMOSZoneA.Range("ZoneATrnCstStart", _ 
wsMOSZoneA.Range("ZoneATrnCstStart").End(xlDown)) 

Set rTrainMOSAList = wsTrainCostZoneA.Range("TrnMOSZoneAListStart", _ 
wsTrainCostZoneA.Range("TrnMOSZoneAListStart") _ 

.End(xlDown)) 

Set rTrainMOSACostList = wsTrainCostZoneA.Range("TrnMOSZoneACostStart", _ 
wsTrainCostZoneA.Range("TrnMOSZoneACostStart") _ 
.End(xlDown)) 

Call CheckMOS (rZoneAMOS, rZoneAMOSTrnVal, rTrainMOSAList, _ 
rTrainMOSACostList) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are MOS(s) in Zone A not listed in the Training" _ 

& " Cost Data-Zone A worksheet." & vbCrLf & " See Data -" _ 

& "Zone A Workseet.") 
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End If 

'Check & Fill in Zone B Ranges 

Set rZoneBMOS = wsMOSZoneB.Range("ZoneBMOSStart", wsMOSZoneB.Range _ 
("ZoneBMOSStart").End(xlDown)) 

Set rZoneBMOSTmVal = wsMOSZoneB .Range("ZoneBTmCstStart", _ 
wsMOSZoneB .Range("ZoneBTmCstStart").End(xlDown)) 

Set rTrainMOSBCList = wsTrainCostZoneBC.Range("TmMOSZoneBCListStart", _ 
wsTrainCostZoneBC.Range("TrnMOSZoneBCListStart") _ 
.End(xlDown)) 

Set rTrainMOSBCCostList = wsTrainCostZoneBC.Range _ 
("TrnMOSZoneBCCostStart", _ 
wsTrainCostZoneBC.Range _ 

("TrnMOSZoneBCCostStart").End(xlDown)) 

Call CheckMOS(rZoneBMOS, rZoneBMOSTmVal, rTrainMOSBCList, _ 
rTrainMOSBCCostList) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are MOS(s) in Zone B not listed in the Training" _ 

& " Cost Data-Zone BC worksheet." & vbCrLf & " See Data -" _ 

& "Zone B Workseet.") 

End If 

'Check & Fill in Zone C Ranges 

Set rZoneCMOS = wsMOSZoneC.Range("ZoneCMOSStart", wsMOSZoneC.Range _ 
("ZoneCMOSStart").End(xlDown)) 

Set rZoneCMOSTmVal = wsMOSZoneC.Range("ZoneCTmCstStart", _ 
wsMOSZoneC. Range("ZoneCTmCstStart").End(xlDown)) 

Call CheckMOS(rZoneCMOS, rZoneCMOSTmVal, rTrainMOSBCList, _ 
rTrainMOSBCCostList) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are MOS(s) in Zone C not listed in the Training" _ 

& " Cost Data-Zone BC worksheet." & vbCrLf & " See Data -" _ 

& "Zone C Workseet.") 

End If 
End Sub 

’CheckMOS() is a subroutine that will verify that the MOS is listed in the 
'Training Cost Worksheet. If the MOS is listed the training cost will be listed 
'in the training cost column, if not an error will be written. 

f 

'@param - rMOSToCheck is the range ofMOSs in the worksheet that need to be verified 

- rTrnValToFill is the range where training costs will be filled in the 
worksheet 

- rMOSList is the range ofMOSs that are listed in the Traing Cost Worksheet 

' - rTrnCost is the range of training costs for MOSs listed in the Traing Cost 

Worksheet 
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Sub CheckMOS(rMOSToCheck As Range, rTrnValToFill As Range, _ 
rMOSList As Range, rTmCost As Range) 

Dim INumMOS As Long 
Dim INumMOSToCheck As Long 
INumErrors = 0 
INumNodes = 0 

INumMOS = rMOSList.Rows.Count 
INumMOSToCheck = rMOSToCheck.Rows.Count 
IMaxIndex = 3 * INumMOS 
ReDim llndex(l To IMaxIndex) As Long 
ReDim sNodeNames(0 To 2 * INumMOS) As String 
Dim i As Integer 
For i = 1 To INumMOS 
AddNode rMOSList(i) 

Next i 

rTmValToFill.Clear 

For i = 1 To INumMOSToCheck 

If DoesContain(rMOSToCheck(i)) Then 
rTmValToFill(i) = rTrnCost(i) 

Else 

rTmValToFill(i) = "Error, MOS not LISTED" 

INumErrors = INumErrors + 1 
End If 
Next i 

Call FormatTwoDecimal(rTrnValToFill) 

End Sub 

ReenlistmentRateCheckByOccFld() is a subroutine that will verify that the OccField 
'is listed in the Reenlistment Rate Worksheet for the correct zone. This 
'subroutine will call ReenlistmentrateCheckMaster and make it verify by OccField. 

t 

'@param - None 

f 

Sub ReenlistmentRateCheckB yOccFld() 

Call ReenlistmentRateCheckMaster(True) 

End Sub 

'ReenlistmentRateCheckByMOS() is a subroutine that will verify that the MOS is listed 
'in the Reenlistment Rate Worksheet for the correct zone. This subroutine will call 
'ReenlistmentrateCheckMaster and make it verify by MOS. 

f 

'@param - None 

f 

Sub ReenlistmentRateCheckByMOS() 

Call ReenlistmentRateCheckMaster(False) 
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End Sub 


'ReenlislmenlRaleCheckMasleii) is a subroutine that will verify that the MOS/OccField 
'is listed in the Reenlistment Rate Worksheet for the correct zone. This subroutine 
'gathers the data to use and calls the subroutine CheckRRR() to verify that the 
'MOS/OccField is listed. If the MOS/OccField is listed the reenlistment rates and its 
'corresponding multiplier will be listed in the reenlistment rate column, if not an 
'error will be written and the user will be notified. 

t 

'@parcim - CheckByOccField is a boolean variable that will be True if the reenlistment 
rates need to be verified by OccFielcl or false if the rates need to be 
verified by MOS 

t 

Sub ReenlistmentRateCheckMaster(bCheckByOccField As Boolean) 

Dim lLength As Long 
Dim lWidth As Long 
'Data Sheet Ranges 
Dim rZoneAMOS As Range 
Dim rZoneAOccField As Range 
Dim rZoneAEASPop As Range 
Dim rZoneAB Space As Range 
Dim rZoneARRR As Range 
Dim rZoneBMOS As Range 
Dim rZoneBOccField As Range 
Dim rZoneBEASPop As Range 
Dim rZoneBBSpace As Range 
Dim rZoneBRRR As Range 
Dim rZoneCMOS As Range 
Dim rZoneCOccField As Range 
Dim rZoneCEASPop As Range 
Dim rZoneCBSpace As Range 
Dim rZoneCRRR As Range 
'Reenlistment Rates Ranges 
Dim rRRRZoneAOccField As Range 
Dim rRRRZoneBOccField As Range 
Dim rRRRZoneCOccField As Range 
Dim rRRRZoneA As Range 
Dim rRRRZoneB As Range 
Dim rRRRZoneC As Range 
'Check Zone A by OccFielcl 

Set rZoneAMOS = wsMOSZoneA.Range("ZoneAMOSStart", wsMOSZoneA.Range 
("ZoneAMOSStart").End(xlDown)) 
lLength = rZoneAMOS.Rows.Count 

Set rZoneAOccField = wsMOSZoneA.Range("ZoneAOccFieldStart", _ 
wsMOSZoneA.Range("ZoneAOccFieldStart") _ 
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.Offset(lLength, 0)) 

Set rZoneAEASPop = wsMOSZoneA.Range("ZoneAEASPopStart", _ 
wsMOSZoneA.Range("ZoneAEASPopStart") _ 

.End(xlDown)) 

Set rZoneABSpace = wsMOSZoneA.Range("ZoneABSAdjStart", _ 
wsMOSZoneA.Range("ZoneABSAdjStart") _ 

.End(xlDown)) 

Set rRRRZoneAOccField = wsRRRZoneA.Range("RRRZoneAOccFieldStart", _ 
wsRRRZoneA.Range("RRRZoneAOccFieldStart") _ 

.End(xlDown)) 

Set rRRRZoneA = wsRRRZoneA.Range(wsRRRZoneA.Range("RRRZoneAStart") 
.End(xlToRight), wsRRRZoneA.Range("RRRZoneAStart") _ 
.End(xlDown)) 

lWidth = rRRRZoneA.Columns.Count 

Set rZoneARRR = wsMOSZoneA.Range("ZoneARRRStart", wsMOSZoneA.Range 
("ZoneARRRStart").Offset(lLength + 1, lWidth + 1)) 

Call CheckRRR(rZoneAOccField, rZoneAMOS, rZoneAEASPop, rZoneABSpace, _ 
rZoneARRR, rRRRZoneAOccField, rRRRZoneA, lLength, _ 
lWidth, bCheckByOccField) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are OccField/MOS(s) in Zone A not listed in the " _ 

& "Reenlistment Rates-Zone A worksheet." & vbCrLf & " See Data - " _ 

& "Zone A Workseet.") 

End If 

'Check Zone B by OccFielcl 

Set rZoneBMOS = wsMOSZoneB.Range("ZoneBMOSStart", wsMOSZoneB.Range 
("ZoneBMOSStart").End(xlDown)) 
lLength = rZoneBMOS. Rows.Count 

Set rZoneBOccField = wsMOSZoneB.Range("ZoneBOccFieldStart", _ 
wsMOSZoneB.Range("ZoneBOccFieldStart") _ 

.Offset(lLength, 0)) 

Set rZoneBEASPop = wsMOSZoneB.Range("ZoneBEASPopStart", _ 
wsMOSZoneB.Range("ZoneBEASPopStart") __ 

.End(xlDown)) 

Set rZoneBBSpace = wsMOSZoneB.Range("ZoneBBSAdjStart", __ 
wsMOSZoneB.Range("ZoneBBSAdjStart") _ 

.End(xlDown)) 

Set rRRRZoneBOccField = wsRRRZoneB.Range("RRRZoneBOccFieldStart", _ 
wsRRRZoneB.Range("RRRZoneBOccFieldStart") _ 

.End(xlDown)) 

Set rRRRZoneB = wsRRRZoneB.Range(wsRRRZoneB.Range("RRRZoneBStart") _ 
.End(xlToRight), wsRRRZoneB.Range("RRRZoneBStart") _ 
.End(xlDown)) 

lWidth = rRRRZoneB.Columns.Count 

Set rZoneBRRR = wsMOSZoneB.Range("ZoneBRRRStart", wsMOSZoneB.Range _ 
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("ZoneBRRRStart").Offset(lLength + 1, lWidth + 1)) 

Call CheckRRR(rZoneBOccField, rZoneBMOS, rZoneBEASPop, rZoneBBSpace, _ 
rZoneBRRR, rRRRZoneBOccField, rRRRZoneB, lLength, _ 
lWidth, bCheckByOccField) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are OccField/MOS(s) in Zone B not listed in the " _ 

& "Reenlistment Rates-Zone B worksheet." & vbCrLf & " See Data - " 

& "Zone B Workseet.") 

End If 

'Check Zone C by OccField 

Set rZoneCMOS = wsMOSZoneC.Range("ZoneCMOSStart", wsMOSZoneC.Range 
("ZoneCMOSStart").End(xlDown)) 
lLength = rZoneCMOS.Rows.Count 

Set rZoneCOccField = wsMOSZoneC.Range("ZoneCOccFieldStart", _ 
wsMOSZoneC.Range("ZoneCOccFieldStart") _ 

.Offset(lLength, 0)) 

Set rZoneCEASPop = wsMOSZoneC.Range("ZoneCEASPopStart", _ 
wsMOSZoneC.Range("ZoneCEASPopStart") _ 

.End(xlDown)) 

Set rZoneCBSpace = wsMOSZoneC.Range("ZoneCBSAdjStart", _ 
wsMOSZoneC.Range("ZoneCBSAdjStart") _ 

.End(xlDown)) 

Set rRRRZoneCOccField = wsRRRZoneC.Range("RRRZoneCOccFieldStart", _ 
wsRRRZoneC.Range("RRRZoneCOccFieldStart") _ 

.End(xlDown)) 

Set rRRRZoneC = wsRRRZoneC.Range(wsRRRZoneC.Range("RRRZoneCStart") _ 
.End(xlToRight), wsRRRZoneC.Range("RRRZoneCStart") _ 
.End(xlDown)) 

lWidth = rRRRZoneC.Columns.Count 

Set rZoneCRRR = wsMOSZoneC.Range("ZoneCRRRStart", wsMOSZoneC.Range _ 
("ZoneCRRRStart").Offset(lLength + 1, lWidth + 1)) 

Call CheckRRR(rZoneCOccField, rZoneCMOS, rZoneCEASPop, rZoneCBSpace, _ 
rZoneCRRR, rRRRZoneCOccField, rRRRZoneC, lLength, __ 
lWidth, bCheckByOccField) 

If INum E rrors > 0 Then 

MsgBox ("Error: There are OccField/MOS(s) in Zone C not listed in the " _ 

& "Reenlistment Rates-Zone C worksheet." & vbCrLf & " See Data -" _ 

& "Zone C Workseet.") 

End If 
End Sub 

'CheckRRR() is a subroutine that will verify that the MOS/OccFielcl is listed in the 
Reenlistment Rate Worksheet for the correct zone. If the MOS/OccFielcl is listed 
'the reenlistment rates and its corresponding multiplier will be listed in the 
'reenlistment rate column, if not an error will be written. 
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'@param - rOccFieldZoneList is the range where the correct OccField will be listed 
on the Data Sheet 

- rZoneMOSList is the range ofMOS's that need to be verified 

- rEASPop is the range ofEAS Populations corresponding to the MOS that is 
to be verified 

' - rBSpace is the range of Boat Spaces available for the MOS that is to be 

verified 

- rZoneRRR is the range that will be filled with reenlistment rates from 
the Reenlistment Rate Worksheet 

- rRRRZoneList is the range of MOS or OccFielcls that the ZoneMOSList will be 
verified against 

- rRRRZone is the range of re enlistment rates corresponding to the RRRZoneList 

- INumToCheck is the number of MOS/OccFields to check 

- INumMult is the number of SRB multipliers that are being used 

- bCheckByOccField is a boolean variable that will be True if the reenlistment 
rates need to be verified by OccField or false if the rates need to be 
verified by MOS 

t 

Sub CheckRRR(rOccFieldZoneList As Range, rZoneMOSList As Range, __ 
rEASPop As Range, rBSpace As Range, rZoneRRR As Range, _ 
rRRRZoneList As Range, rRRRZone As Range, INumToCheck As Long, _ 
INumMult As Long, bCheckByOccLield As Boolean) 

Dim INumRRR As Long 
Dim lPoint As Long 
Dim sOccLield As String 
Dim bContain As Boolean 
rOccLieldZoneList. Clear 
rZoneRRR. Clear 
INumErrors = 0 
INumNodes = 0 

INumRRR = rRRRZoneList.Rows.Count 

IMaxIndex = 3 * INumRRR 

ReDim llndex(l To IMaxIndex) As Long 

ReDim sNodeNames(0 To 2 * INumRRR) As String 

Dim i As Integer 

Dim j As Integer 

Dim lRRRPointer() As Long 

ReDim lRRRPointer(l To IMaxIndex) 

Lor i = 1 To INumRRR 
AddNode rRRRZoneList(i) 

lRRRPointer(LookupNode(rRRRZoneList(i))) = i + 1 
Next i 

Formating here so the 0's on the MOSs will show up 
Call LormatOccLieldText(rOccLieldZoneList) 
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'Prints the multiplier Numbers 
For i = 1 To INumMult 

rZoneRRR(l, i) = rRRRZone(l, i) 

Next i 

For i = 1 To INumToCheck 

sOccField = Mid(rZoneMOSList(i), 1, 2) 

If bCheckByOccField Then 

bContain = DoesContain(sOccField) 

If bContain Then 

lPoint = lRRRPointer(LookupNode(sOccField)) 

End If 
Else 

bContain = DoesContain(rZoneMOSList(i)) 

If bContain Then 

lPoint = lRRRPointer(LookupNode(rZoneMOSList(i))) 

End If 
End If 

rOccFieldZoneList(i) = sOccField 
If bContain Then 

For j = 1 To INumMult 

rZoneRRR(i + 1, j) = Application.WorksheetFunction.Round _ 
(rRRRZone(lPoint, j), 4) 

Next j 
Else 

If rEASPop(i) = 0 Then 

rZoneRRR(i + 1, 1) = "Error, MOS not LISTED and NOT needed due " 
& "to Zero EAS Population" 

Elself rBSpace(i) = 0 Then 

rZoneRRR(i + 1, 1) = "Error, MOS not LISTED and NOT needed due " 
& "to Zero Boat Spaces" 

Else 

rZoneRRR(i +1,1) = "Error, MOS not LISTED and NEEDED" 
INumErrors = INumErrors + 1 
End If 
End If 
Next i 

Call FormatTwoDecimal(rZoneRRR) 

End Sub 

'File: HashCode Module in VBA 
'Date Created: March 2007 
'Last Updated: May 2007 

f 

'Checks is a VBA Module that contains a list of subroutines that will perform 
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'checks on the data set to verify that it is complete in order to run the 
'optimization software. 

t 

'@ Author Matthew Caryle 
' Updated by Kent Robbins, Jr. 

t 

Option Explicit 

'HashName is a function that will take a string variable and convert it to a 
'numerical value for use in a hash code. 

f 

'@ par am - sName is a String Variable that will be converted to a numerical value 
'@ return - HashName is a Long variable that is converted from sName 

f 

Function HashName(sName As String) As Long 

'Converts an input string sName into a long integer between 1 and IMaxIndex 
Dim lValue As Long 
Dim i As Integer 
lValue = 0 

For i = 1 To Len(sName) 

lValue = (10 * lValue + Asc(Mid(sName, i, 1))) Mod IMaxIndex 
Next i 

HashName = lValue + 1 
End Function 

AddNode is a subroutine that will add a variable to the hash code. 

f 

'@ par am - sName is a string variable that will be added to the hash code. 

t 

Sub AddNode(sName As String) 

Dim hv As Long 

hv = HashName(sName) 

Do While llndex(hv) <> 0 And sNodeNames(lIndex(hv)) <> sName 
hv = hv + 1 

If hv > IMaxIndex Then 
hv = 1 
End If 
Loop 

If llndex(hv) = 0 Then 

INumNodes = INumNodes + 1 
llndex(hv) = INumNodes 
sNodeNames(lNumNodes) = sName 
End If 
End Sub 

'LoopupNode is a function that will return the node number of the parameter. 
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'@param - sName is a string variable used to look up its associated node number. 
'@ return - LoopUpNode 


Function LookupNode(sName As String) As Long 
Dim hv As Long 
hv = HashName(sName) 

Do While llndex(hv) <> 0 And sNodeNames(Hndex(hv)) <> sName 
hv = hv + 1 

If hv > IMaxIndex Then 
hv = 1 
End If 
Loop 

LookupNode = llndex(hv) 

End Function 

'DoesContain is a function that will compute the correct number of reenlistments for 
'Zone A that are needed when BNA school seats are used. 

f 

'@param - sName is a string variable that needs to be check to see if it is in the 
1 hash code. 

'@ return - DoesContain is a boolean variable. True if sName is in the hash code, 
otherwise false. 

f 

Function DoesContain(sName As String) As Boolean 
Dim hv As Long 
hv = HashName(sName) 

Do While llndex(hv) <> 0 And sNodeNames(Hndex(hv)) <> sName 
hv = hv + 1 

If hv > IMaxIndex Then 
hv = 1 
End If 
Loop 

If llndex(hv) = 0 Then 
DoesContain = False 
Else 

DoesContain = True 
End If 

End Function 

'File: Reduction Module in VBA 
'Date Created: 17 June 2007 
'Last Updated: 25 July 2007 

t 

'Reduction is a VBA Modula that will combine all MOSfrom Zone A, B, and C onto a 
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'single worksheet containing the necessary data for the GAMS and VBA solver to find 
'the Optimum SRB multiplier. Reduction will only transfer MOSs that will not achieve 
'the reenlistment goal without the use of a multiplier greater than 0. Reduction 
'will also check the cost of each multiplier and will assign it a 0.0 reenlistment rate 
'if that multiplier violates the maximum bonus contstraint. 

f 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 

'ReduceSizeMaster is a subroutine that will gather the data neccessary and call the 
function Reduce to execute the reduction and transfer of the MOSsin Zones A, B, and 
'C that will not achieve the reenlistment goal with a multiplier ofO. Results are 
'transferecl to the ReducedFull Problem Worksheet. 

t 

'@param - none 

f 

Sub ReduceSizeMaster() 

'Zone Data Ranges 
Dim rZoneAOccField As Range 
Dim rZoneAMOSField As Range 
Dim rZoneAZoneField As Range 
Dim rZoneAMOSPopField As Range 
Dim rZoneAEASPopField As Range 
Dim rZoneAPreMultField As Range 
Dim rZoncAWcightFicld As Range 
Dim rZoneAReupNeedField As Range 
Dim rZoneATrainCostField As Range 
Dim rZoneARRR As Range 
Dim dZoneAAvgPay As Double 
Dim rZoneBOccField As Range 
Dim rZoneBMOSField As Range 
Dim rZoncBZoncFicld As Range 
Dim rZoneBMOSPopField As Range 
Dim rZoneBEASPopField As Range 
Dim rZoneBPreMultField As Range 
Dim rZoneBWeightField As Range 
Dim rZoneBReupNeedField As Range 
Dim rZoneBTrainCostField As Range 
Dim rZoneBRRR As Range 
Dim dZoneBAvgPay As Double 
Dim rZoneCOccField As Range 
Dim rZoneCMOSField As Range 
Dim rZoneCZoneField As Range 
Dim rZoneCMOSPopField As Range 
Dim rZoneCEASPopField As Range 
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Dim rZoneCPreMultField As Range 
Dim rZoneCWeightField As Range 
Dim rZoneCReupNeedField As Range 
Dim rZoneCTrainCostField As Range 
Dim rZoneCRRR As Range 
Dim dZoneCAvgPay As Double 
Dim dMaxIndivBonus As Double 
'OutPutZone Variables 
Dim iMaxOutputLength As Integer 
Dim iMaxOutputWidth As Integer 
Dim rOutput As Range 
Dim rAvgPays As Range 
Dim i As Integer 
Dim iLinePosition As Integer 

dMaxIndivBonus = wsReducedProblem.Range("MaxIndBonus") 

Set rZoneAOccField = wsMOSZoneA.Range("ZoneAOccFieldStart", _ 
wsMOSZoneA.Range("ZoneAOccFieldStart") _ 

.End(xlDown)) 

Set rZoneAMOSField = wsMOSZoneA.Range("ZoneAMOSStart", _ 
wsMOSZoneA.Range("ZoneAMOSStart") _ 

.End(xlDown)) 

Set rZoneAZoneField = wsMOSZoneA.Range("ZoneASRBZoneStart", _ 
wsMOSZoneA.Range("ZoneASRBZoneStart") _ 

.End(xlDown)) 

Set rZoneAMOSPopField = wsMOSZoneA.Range("ZoneAMOSPopStart", __ 
wsMOSZoneA.Range("ZoneAMOSPopStart") _ 

.End(xlDown)) 

Set rZoneAEASPopField = wsMOSZoneA.Range("ZoneAEASPopStart", _ 
wsMOSZoneA.Range("ZoneAEASPopStart") _ 

.End(xlDown)) 

Set rZoneAPreMultField = wsMOSZoneA.Range("ZoneAPresetStart", _ 
wsMOSZoneA.Range("ZoneAPresetStart") _ 

.End(xlDown)) 

Set rZoneAWeightField = wsMOSZoneA.Range("ZoneAWgtFctStart", _ 
wsMOSZoneA.Range("ZoneAWgtFctStart") _ 

.End(xlDown)) 

Set rZoneAReupNeedField = wsMOSZoneA.Range("ZoneABSAdjStart", _ 
wsMOSZoneA.Range("ZoneABSAdjStart") _ 

.End(xlDown)) 

Set rZoneATrainCostField = wsMOSZoneA.Range("ZoneATmCstStart", _ 
wsMOSZoneA.Range("ZoneATmCstStart") _ 

.End(xlDown)) 

Set rZoneARRR = wsMOSZoneA.Range(wsMOSZoneA.Range("ZoneARRRStart") 
.End(xlToRight), wsMOSZoneA.Range("ZoneARRRStart") _ 
.End(xlDown)) 
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dZoneAAvgPay = wsMOSZoneA.Range("ZoneAAvgPay") 

Set rZoneBOccField = wsMOSZoneB.Range("ZoneBOccFieldStart", _ 
wsMOSZoneB.Range("ZoneBOccFieldStart") _ 

.End(xlDown)) 

Set rZoneBMOSField = wsMOSZoneB.Range("ZoneBMOSStart", _ 
wsMOSZoneB.Range("ZoneBMOSStart") _ 

.End(xlDown)) 

Set rZoneBZoneField = wsMOSZoneB.Range("ZoneBSRBZoneStart", _ 
wsMOSZoneB .Range("ZoneBSRBZoneStart") _ 

.End(xlDown)) 

Set rZoneBMOSPopField = wsMOSZoneB .Range("ZoneBMOSPopStart", _ 
wsMOSZoneB .Range("ZoneBMOSPopStart") _ 

.End(xlDown)) 

Set rZoneBEASPopField = wsMOSZoneB.Range("ZoneBEASPopStart", _ 
wsMOSZoneB .Range("ZoneBEASPopStart") _ 

.End(xlDown)) 

Set rZoneBPreMultField = wsMOSZoneB.Range("ZoneBPresetStart", __ 
wsMOSZoneB.Range("ZoneBPresetStart") _ 

.End(xlDown)) 

Set rZoneBWeightField = wsMOSZoneB.Range("ZoneBWgtFctStart", _ 
wsMOSZoneB .Range("ZoneBWgtFctStart") _ 

.End(xlDown)) 

Set rZoneBReupNeedField = wsMOSZoneB .Range("ZoneBBSAdjStart", __ 
wsMOSZoneB .Range("ZoneBBSAdjStart") _ 

.End(xlDown)) 

Set rZoneBTrainCostField = wsMOSZoneB.Range("ZoneBTrnCstStart", _ 
wsMOSZoneB.Range("ZoneBTrnCstStart") _ 

.End(xlDown)) 

Set rZoneBRRR = wsMOSZoneB.Range(wsMOSZoneB.Range("ZoneBRRRStart") 
.End(xlToRight), wsMOSZoneB.Range("ZoneBRRRStart") _ 
.End(xlDown)) 

dZoneBAvgPay = wsMOSZoneB.Range("ZoneBAvgPay") 

Set rZoneCOccField = wsMOSZoneC.Range("ZoneCOccFieldStart", _ 
wsMOSZoneC.Range("ZoneCOccFieldStart") _ 

.End(xlDown)) 

Set rZoneCMOSField = wsMOSZoneC.Range("ZoneCMOSStart", _ 
wsMOSZoneC.Range("ZoneCMOSStart") _ 

.End(xlDown)) 

Set rZoneCZoneField = wsMOSZoneC.Range("ZoneCSRBZoneStart", _ 
wsMOSZoneC.Range("ZoneCSRBZoneStart") _ 

.End(xlDown)) 

Set rZoneCMOSPopField = wsMOSZoneC.Range("ZoneCMOSPopStart", _ 
wsMOSZoneC.Range("ZoneCMOSPopStart") _ 

.End(xlDown)) 

Set rZoneCEASPopField = wsMOSZoneC.Range("ZoneCEASPopStart", _ 
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wsMOSZoneC.Range("ZoneCEASPopStart") _ 

.End(xlDown)) 

Set rZoneCPreMultField = wsMOSZoneC.Range("ZoneCPresetStart", _ 
wsMOSZoneC.Range("ZoneCPresetStart") _ 

.End(xlDown)) 

Set rZoneCWeightField = wsMOSZoneC.Range("ZoneCWgtFctStart", _ 
wsMOSZoneC.Range("ZoneCWgtFctStart") _ 

.End(xlDown)) 

Set rZoneCReupNeedField = wsMOSZoneC.Range("ZoneCBSAdjStart", __ 
wsMOSZoneC.Range("ZoneCBSAdjStart") _ 

.End(xlDown)) 

Set rZoneCTrainCostField = wsMOSZoneC.Range("ZoneCTrnCstStart", _ 
wsMOSZoneC.Range("ZoneCTrnCstStart") _ 

.End(xlDown)) 

Set rZoneCRRR = wsMOSZoneC.Range(wsMOSZoneC.Range("ZoneCRRRStart") 
.End(xlToRight), wsMOSZoneC.Range("ZoneCRRRStart") _ 
.End(xlDown)) 

dZoneCAvgPay = wsMOSZoneC.Range("ZoneCAvgPay") 

iMaxOutputLength = rZoneAOccField.Rows.Count + rZoneBOccField.Rows.Count 
+ rZoneCOccField.Rows.Count 
iMaxOutputWidth = 9 + rZoneARRR.Columns.Count 
iLinePosition = 2 

Set rOutput = wsReducedProblem.Range("RFPOutputStart", __ 
wsReducedProblem.Range("RFPOutputStart") _ 

.Offset(iMaxOutputLength + 1, iMaxOutputWidth)) 

Set rAvgPays = wsReducedProblem.Range("AvgPayA", wsReducedProblem.Range 
("AvgPay A") .Offset(3, 0)) 
rOutput.Clear 

Call FormatOccFieldText(rOutput) 
rOutput.HorizontalAlignment = xlCenter 
rAvgPays.Clear 

rAvgPays.HorizontalAlignment = xlCenter 

rOutput(l, 1) = "Occ Field" 

rOutput(l, 2) = "MOS" 

rOutput(l, 3) = "Zone" 

rOutput(l, 4) = "MOSPop" 

rOutput(l, 5) = "EASPop" 

rOutput(l, 6) = "PresetMult" 

rOutput(l, 7) = "WgtFactor" 

rOutput(l, 8) = "B/SPACE" 

rOutput(l, 9) = "TrnCosts" 

For i = 1 To rZoneARRR.Columns.Count 
rOutput(l, 9 + i) = rZoneARRR(l, i) 

Next i 

rAvgPays(l) = dZone A AvgPay 
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rAvgPays(2) = dZoneBAvgPay 
rAvgPays(3) = dZoneCAvgPay 

iLinePosition = Reduce(rZoneAOccField, rZoneAMOSField, rZoneAZoneField, _ 
rZoneAMOSPopField, rZoneAEASPopField, rZoneAPreMultField, _ 
rZoneAWeightField, rZoneAReupNeedField, rZoneATrainCostField, _ 
rZoneARRR, rOutput, dZoneAAvgPay, dMaxIndivBonus, iLinePosition) 
iLinePosition = Reduce(rZoneBOccLield, rZoneBMOSLield, rZoneBZoneLield, _ 
rZoneBMOSPopLield, rZoneBEASPopLield, rZoneBPreMultLield, __ 
rZoneBWeightLield, rZoneBReupNeedLield, rZoneBTrainCostLield, _ 
rZoneBRRR, rOutput, dZoneBAvgPay, dMaxIndivBonus, iLinePosition) 
iLinePosition = Reduce(rZoneCOccLield, rZoneCMOSLield, rZoneCZoneLield, _ 
rZoneCMOSPopLield, rZoneCEASPopLield, rZoneCPreMultLield, _ 
rZoneCWeightLield, rZoneCReupNeedLield, rZoneCTrainCostLield, _ 
rZoneCRRR, rOutput, dZoneCAvgPay, dMaxIndivBonus, iLinePosition) 

End Sub 

'Reduce is a function that will execute the reduction and transfer of those MOSs in 
'a Zone that will need a multiplier greater than 0 to possible achieve the reenlistment 
'goal. Reduce will also check the cost of using each multiplier and if the cost is 
'greater than the maximum allowable bonus the reenlistment rate for that multiplier will 
'be set to 0 so it will not be used. 

f 

'@ par am - rOccField is the range of OccFields of the Zones MOSs to be reduced. 

- rMos is the range of MOSs in the Zone that is to reduced. 

- rZone is a range variable and is the Zone that to be reduced for each MOS. 

- rMosPop is a range variable of MOS populations for each MOS. 

- rEASPop is a range variable ofEAS populations for each MOS. 

- rPreset is a range variable of preset multiplier value for each MOS. 

- rWeight is a range variable of exogenous weighting factor for each MOS. 

- rBSpace is a range variable of needed reenlistment numbers for each MOS. 

- rTrnCst is a range variable of training costs for each MOS. 

- rRRR is a range variable of reenlistment response rates for each multiplier. 

- rOutput is a range where the MOS and data will be transfered to. 

- dAvgPay is a double variable of the average pay of the Zone being reduced. 

- dMaxBonus is a double variable of the maximum allowable bonus. 

- iLine is an integer variable and a place holder for the location in the 
output for the next Zone to start being transfered to. 

f 

'@ return - Reduce is an integer variable and is the location marker for the output 
for the next Zone that is to be reduced. 

f 

Lunction Reduce(rOccLield As Range, rMos As Range, rZone As Range, _ 
rMosPop As Range, rEASPop As Range, rPreset As Range, _ 
rWeight As Range, rBSpace As Range, rTrnCst As Range, _ 
rRRR As Range, rOutput As Range, dAvgPay As Double, _ 


44 




dMaxBonus As Double, iLine As Integer) As Integer 
Dim dGap As Double 
Dim i As Integer 
Dim j As Integer 
Dim bContinue As Boolean 
dGap = rRRR(l, 2) - rRRR(l, 1) 

For i = 1 To rOccField.Rows.Count 
If rEASPop(i) = 0 Then 
bContinue = False 
Elself rBSpace(i) = 0 Then 
bContinue = False 
Else 

bContinue = True 
End If 

If bContinue Then 
If rPreset(i) > -1 Then 

rOutput(iLine, 1) = rOccField(i) 
rOutput(iLine, 2) = rMos(i) 
rOutput(iLine, 3) = rZone(i) 
rOutput(iLine, 4) = rMosPop(i) 
rOutput(iLine, 5) = rEASPop(i) 
rOutput(iLine, 6) = rPreset(i) 
rOutput(iLine, 7) = rWeight(i) 
rOutput(iLine, 8) = rBSpace(i) 
rOutput(iLine, 9) = rTrnCst(i) 

For j = 1 To rRRR.Columns.Count 
rOutput(iLine, 9 + j) = rRRR(i + 1, j) 

Next j 

iLine = iLine + 1 

Elself (rEASPop(i) * rRRR(i + 1, 1) < rBSpace(i)) Then 
rOutput(iLine, 1) = rOccField(i) 
rOutput(iLine, 2) = rMos(i) 
rOutput(iLine, 3) = rZone(i) 
rOutput(iLine, 4) = rMosPop(i) 
rOutput(iLine, 5) = rEASPop(i) 
rOutput(iLine, 6) = rPreset(i) 
rOutput(iLine, 7) = rWeight(i) 
rOutput(iLine, 8) = rBSpace(i) 
rOutput(iLine, 9) = rTrnCst(i) 

Forj = 1 To rRRR.Columns.Count 

If (dGap * (j - 1) * dAvgPay * 4 > dMaxBonus) Then 
rOutput(iLine, 9 + j) = 0 
Else 

rOutput(iLine, 9 + j) = rRRR(i + 1, j) 

End If 
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Next j 

iLine = iLine + 1 
End If 
End If 
Next i 

Reduce = iLine 
End Function 

'File: SolveVBA Module in VBA 
'Date Created: 17 June 2007 
'Last Updated: 31 August 2007 

f 

'SolveVBA Module is a VBA module that will solve the Lagrangian relaxation problem of 
finging the optimum SRB multipliers to assign to each MOSZ based on the weighted 
'penalty function occured with each multiplier and subject to certain buget 
'limitations. 

f 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 
Public dlnfinity As Double 
Public dBudget As Double 
Public dTest As Double 
Public iNumMOS As Integer 
Public iNumMult As Integer 
Public dGap As Double 
Public dHLambda As Double 
Public dZL As Double 
Public dZU As Double 
Public dZH As Double 
Public dSumCost As Double 
Public dFinalCost As Double 
Public dTotalDeviation As Double 

'SolveEXCEL is the master subroutine that run the VBA Model. This subroutine will 
'gather the pertinent data that is need in later subroutines to find the optimum 
'SRB multipliers. 

f 

'@param - none 

t 

Sub SolveEXCEL() 

Dim dStartTime As Double 
dStartTime = Timer 
Dim dFinishTime As Double 
'Worksheet variables 
Dim rMOSData As Range 


46 




Dim rBudgetlnfo As Range 
Dim dCost() As Double 
Dim dDeviation() As Double 
Dim i As Integer 
Dim dSoln() As Double 
Dim dPreset() As Double 
Dim dOptimality As Double 
dTotalDeviation = 0# 

'Read in data from worksheets 

Set rMOSData = wsReducedProblem.Range(wsReducedProblem.Range _ 
("RFPOutputStart").End(xlToRight), wsReducedProblem.Range _ 
("RFPOutputStart").End(xlDown)) 

Set rBudgetlnfo = wsReducedProblem.Range("OverageWeight", _ 

wsReducedProblem.Range("OverageWeight").End(xlDown)) 
dBudget = rBudgetInfo(2) 
dlnfinity = 1 * (10 A 22) 
iNumMOS = rMOSData.Rows.Count 
iNumMult = rMOSData.Columns.Count - 9 
dGap = rMOSData(l, 11) - rMOSData(l, 10) 

ReDim dCost(l To (iNumMOS - 1), 1 To iNumMult) As Double 
ReDim dDeviation(l To (iNumMOS - 1), 1 To iNumMult) As Double 
ReDim dPreset(l To (iNumMOS - 1), 1 To 2) As Double 
ReDim dSoln(l To (iNumMOS - 1)) As Double 
For i = 2 To iNumMOS 

If rMOSData(i, 6) >= 0 Then 

dPreset(i - 1, 1) = rMOSData(i, 6) 
dPreset(i - 1,2) = rMOSData(i, 6) / dGap + 1 
Else 

dPreset(i - 1, 1) = -1 
dPreset(i - 1, 2) = -1 
End If 
Next i 

Call FillCostandDeviation(rMOSData, rBudgetlnfo, dCost(), dDeviation()) 

Call Bound(dSoln(), dPreset(), dCost(), dDeviation()) 

Call Heuris(dDeviation(), dCost(), dSoln(), dPreset()) 

Call WriteResults(rMOSData, rBudgetlnfo, dSoln()) 

dFinishTime = Timer 

dOptimality = (dZH - dZL) / dZL * 100 

MsgBox ("VBA Solve has finished." & vbCrLf & "Solution is within " & Round 
(dOptimality, 6) & "% of optimality" & vbCrLf & "Run time was " & _ 
Round((dFinishTime - dStartTime), 3) & " seconds" & vbCrLf & "Object". 
& "Function is " & Round(dTotalDeviation, 4) & vbCrLf & "Total Cost" _ 
& "is $" & Round(dFinalCost * 1000, 2)) 

End Sub 

'FillCostandDeviation is a subroutine that will fill the Cost and Deviation arrays 
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for each MOSZ ancl multiplier combination. 

f 

'@param - rMOSData is a range variable that contains data for all the MOSZs. 

- rBudgetlnfo is a range variable that contains the budget and pay 
information needed for the SRB Optimization problem. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dDeviation is an array of doubles that will contain the deviation for 
setting each MOSZ with all available multipliers. 

f 

Sub FillCostandDeviation(rMOSData As Range, rBudgetlnfo As Range, _ 
dCost() As Double, dDeviation() As Double) 

Dim i As Integer 
Dim dScaleFactor As Double 
Dim dWeight As Double 
Dim dEligible As Double 
Dim dNeededNum As Double 
Dim dTotalNum As Double 
Dim dTrainCost As Double 
Dim dAvgPay As Double 
Dim iContractLength As Integer 
Dim dMaxTrainCost As Double 
dScaleFactor = 1000# 
dMaxTrainCost = 0# 
dBudget = dBudget / dScaleFactor 
dTest = dlnfinity / 1.1 
dHLambda = -dlnfinity 
'Set max Training Cost 
For i = 2 To iNumMOS 

dMaxTrainCost = Application.WorksheetFunction.Max(dMaxTrainCost, _ 
rMOSData(i, 9)) 

Next i 

For i = 2 To iNumMOS 

Dim dTotalWeight As Double 

Dim j As Integer 

Dim dDev As Double 

Dim dHL As Double 

dTotalNum = rMOSData(i, 4) 

dEligible = rMOSData(i, 5) 

dWeight = rMOSData(i, 7) 

dNeededNum = rMOSData(i, 8) 

dTrainCost = rMOSData(i, 9) 

dAvgPay = rBudgetInfo(rMOSData(i, 3) + 4) 

iContractLength = rBudgetInfo(3) 

If (dTrainCost < 1#) Then 
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dTrainCost = 50# 

End If 

If (dTotalNum > 0#) Then 

dTotalWeight = dWeight * (dTrainCost / dMaxTrainCost) / dTotalNum 
Else 

dTotalWeight = dWeight * (dTrainCost / dMaxTrainCost) 

End If 

If (dEligible < 1#) Then 
dEligible =1# 

End If 

For j = 1 To iNumMult 

dDev = dNeededNum - dEligible * rMOSData(i, 9 + j) 

If (dDev < 0#) Then 

dDeviation(i - 1, j) = dTotalWeight * rBudgetlnfo(l) * (dDev A 2) 

Else 

dDeviation(i - 1, j) = dTotalWeight * (dDev A 2) 

End If 

If (rMOSData(i, 9 + j) = 0) Then 
dCost(i - 1, j) = dlnfinity 
Else 

dCost(i - 1, j) = ((dEligible * rMOSData(i, 9 + j)) * dAvgPay * _ 
iContractLength * rMOSData(l, 9 + j)) / _ 
dScaleFactor 

End If 

If (j > 1) Then 

dHL = (dDeviation(i - 1, 1) - dDeviation(i - 1, j)) / dCost(i - 1, j) 

If (dHL > dHLambda) Then 
dHLambda = dHL 
End If 
End If 
Next j 
Next i 
End Sub 

'Bound is a subroutine that will determine the upper and lower bounds of the 
'Lagrangian relaxation and call two other subroutines to find feasible solution 
'sets with given Lagrangian multipliers and find the best feasible solution set 
'that has the lowest object function. 

'@param - dSoln is an array of doubles that contains the multipliers solution set. 

- dPreset is an array of doubles the list the preset multipliers if chosen. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dDeviation is an array of doubles that will contain the deviation for 
settingeach MOSZ with all available multipliers. 
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Sub Bound(dSoln() As Double, dPreset() As Double, dCost() As Double,. 
dDeviation() As Double) 

Dim dEps As Double 
Dim dEndL As Double 
Dim dEndR As Double 
Dim dZUBest As Double 
Dim dAmbda As Double 
Dim dFLambda As Double 
Dim i As Integer 
dEps = dHLambda / 100000 
dEndL = 0# 

dEndR = 0.01 * dHLambda 
dZUBest = dlnfinity 
i = 0 
Do 

dAmbda = (dEndL + dEndR) / 2# 

Call MinFcn(dAmbda, dPreset(), dCost(), dDeviation()) 

If (dSumCost <= dBudget) Then 
dEndR = dAmbda 
If (dZU <= dZUBest) Then 
dZUBest = dZU 
dFLambda = dAmbda 
End If 
Exit Do 
End If 

dEndL = dAmbda 
dEndR = 10.01 * dEndR 
i = i + 1 

Loop Until (i = 2) 

Do 

dAmbda = (dEndL + dEndR) / 2# 

Call MinFcn(dAmbda, dPreset(), dCost(), dDeviation()) 

If (dSumCost <= dBudget) Then 
dEndR = dAmbda 
If (dZU <= dZUBest) Then 
dZUBest = dZU 
dFLambda = dAmbda 
End If 
Else 

dEndL = dAmbda 
End If 

Loop Until ((dEndR - dEndL) <= dEps) 

Call MinFeasible(dFLambda, dSoln(), dPreset(), dDeviation(), dCost()) 
End Sub 
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'MinFcn is a subroutine that will determine the minimum objective function for 
'the current Lagrangian multiplier. 

f 

'@param - dAmbda is the Lagrangian multiplier current used. 

- clPreset is an array of doubles the list the preset multipliers if chosen. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dDeviation is an array of doubles that will contain the deviation for 
' setting each MOSZ with cdl available multipliers. 

f 

Sub MinFcn(dAmbda As Double, dPreset() As Double, dCost() As Double, _ 
dDeviation() As Double) 

Dim dCelMin As Double 
Dim dCelTot As Double 
Dim dCelObj As Double 
Dim dCObj As Double 
Dim dZJ As Double 
Dim ilndex As Integer 
Dim i As Integer 
Dim j As Integer 
Dim bContinue As Boolean 
dSumCost = 0# 
dCelTot = 0# 
dZU = 0# 

For i = 1 To (iNumMOS - 1) 
bContinue = True 
Do 

If (dPreset(i, 1) > -1) Then 

dCelMin = dDeviation(i, dPreset(i, 2)) + dAmbda * dCost(i, _ 
dPreset(i, 2)) 

dCObj = dDeviation(i, dPreset(i, 2)) 
ilndex = dPreset(i, 2) 

Exit Do 
End If 

dCelMin = dlnfinity 
ilndex = 0 

For j = 1 To (iNumMult) 

If (dDeviation(i, j) <= dTest) Then 

dZJ = dDeviation(i, j) + dAmbda * dCost(i, j) 
dCelObj = dDeviation(i, j) 

If (dZJ <= dCelMin) Then 
dPreset(i, 2) = j 
dCelMin = dZJ 
dCObj = dCelObj 
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ilndex = j 
End If 
End If 
Next j 

bContinue = False 
Loop Until bContinue = False 
dCelTot = dCelTot + dCelMin 
dZU = dZU + dCObj 
dSumCost = dSumCost + dCost(i, ilndex) 

Next i 

dZL = dCelTot - dAmbda * dBudget 
End Sub 

'MinFeasible is a subroutine that will determine minimum feasible object value 
'solution set with the final Langrangian multiplier. This is the feasible solution 
'set with the lowest objective value. 

t 

'@param - dFLambda is a double variable and is the final chosen Langrangian 
' multiplier. 

- dSoln is an array of doubles that contains the multipliers solution set. 

' - dPreset is an array of doubles the list the preset multipliers if chosen. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dDeviation is an array of doubles that will contain the deviation for 
setting each MOSZ with all available multipliers. 

f 

Sub MinFeasible(dFLambda As Double, dSoln() As Double, dPreset() As Double, 
dDeviation() As Double, dCost() As Double) 

Dim dCelTot As Double 
Dim dCObj As Double 
Dim dCelMin As Double 
Dim dZJ As Double 
Dim dCelObj As Double 
Dim ilndex As Integer 
Dim i As Integer 
Dim j As Integer 
Dim bContinue As Boolean 
dSumCost = 0# 
dCelTot = 0# 
dZU = 0# 

For i = 1 To (iNumMOS - 1) 

Do 

If (dPreset(i, 1) > -1) Then 

dCObj = dDeviation(i, dPreset(i, 2)) 
ilndex = dPreset(i, 2) 
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Exit Do 
End If 

dCelMin = dlnfinity 
ilndex = 0 

For j = 1 To (iNumMult) 

If (dDeviation(i, j) <= dTest) Then 

dZJ = dDeviation(i, j) + dFLambda * dCost(i, j) 
dCelObj = dDeviation(i, j) 

If (dZJ <= dCelMin) Then 
dSoln(i) = j 
dCelMin = dZJ 
dCObj = dCelObj 
ilndex = j 
End If 
End If 
Next j 

bContinue = False 
Loop Until bContinue = False 
dZU = dZU + dCObj 
dSumCost = dSumCost + dCost(i, ilndex) 

Next i 
End Sub 

'Heuris is a subroutine that will improve the solution to the Langranian relaxation 
'of the SRB multiplier problem by using a heuristic to consume any residual budget. 

f 

'@param - dDeviation is an array of doubles that will contain the deviation for 
setting each MOSZ with all available multipliers. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dSoln is an array of doubles that contains the multipliers solution set. 

' - clPreset is an array of doubles the list the preset multipliers if chosen. 

f 

Sub Heuris(dDeviation() As Double, dCost() As Double, dSoln() As Double, _ 
dPreset() As Double) 

Dim dCost 1 As Double 
Dim dCost2 As Double 
Dim i As Integer 

Dim bRoomlnBudget As Boolean 
dCost 1 = dSumCost 
bRoomlnBudget = True 
While bRoomlnBudget = True 
If dCost 1 <= dBudget Then 

dCost2 = NewCostl(dDeviation(), dCost(), dSoln(), dPreset(), dCostl) 

End If 
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If dCost2 = dCostl Then 
bRoomlnBudget = False 
Else 

dCostl = dCost2 
End If 
Wend 
dZH = 0# 
dFinalCost = 0# 

For i = 1 To (iNumMOS - 1) 

dZH = dZH + dDeviation(i, dSoln(i)) 
dFinalCost = dFinalCost + dCost(i, dSoln(i)) 
dTotalDeviation = dTotalDeviation + dDeviation(i, dSoln(i)) 

Next i 
End Sub 

'NewCostl is a function used to calculate the rate of return for each MOSZ based on 
'the rate of improvement to the objective function that can be achieve for each 
'dollar. Newcostl will the pick the MOSZ with the largest ROR and increase the 
'MOSZs solution by one step in the multiplier and calculating the new cost of the 
'solution set. 

f 

'@param - dDeviation is an array of doubles that will contain the deviation for 
setting each MOSZ with all available multipliers. 

- dCost is an array of doubles that will contain the cost for setting each 
MOSZ with all available multipliers. 

- dSoln is an array of doubles that contains the multipliers solution set. 

- dPreset is an array of doubles the list the preset multipliers if chosen. 

- dOldCostl is a double variable that is the current cost of the solution 
set. 

t 

'@ return - NewCostl is a double variable that is the new cost of the new solution 
' set that has been heuristically imporved. 

t 

Function NewCostl(dDeviation() As Double, dCost() As Double, dSoln() As Double, 
dPreset() As Double, dOldCostl As Double) As Double 
Dim dUArray() As Double 
Dim dResid As Double 
Dim dRNum As Double 
Dim dRDen As Double 
Dim dUBest As Double 
Dim i As Integer 
Dim j As Integer 
Dim bContinuelnner As Boolean 
Dim ilndex As Integer 
ReDim dUArray(l To iNumMOS - 1) 
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dResid = dBudget - dOldCostl 
For i = 1 To (iNumMOS - 1) 

Do 'continue do loop 
bContinuelnner = True 
If (dPreset(i, 1) > -1) Then 
dUArray(i) = -dlnfinity 
Exit Do 'exits continue do loop 
End If 

If (dSoln(i) >= iNumMult) Then 
dUArray(i) = -dlnfinity 
Exit Do 

Elself (dDeviation(i, (dSoln(i) + 1)) > dTest) Then 
dUArray(i) = -dlnfinity 
Exit Do 
End If 

dRNum = dDeviation(i, dSoln(i)) - dDeviation(i, (dSoln(i) + 1)) 
dRDen = dCost(i, (dSoln(i) + 1)) - dCost(i, dSoln(i)) 

If ((dRDen > dResid) Or (dRNum <= (0.1 * 10 A -8))) Then 
dUArray(i) = -dlnfinity 
Exit Do 'exits continue do loop 
Else 

dUArray(i) = dRNum / dRDen 
End If 

If (dUArray(i) < 0#) Then 
dUArray(i) = -1 * dUArray(i) 

End If 

Loop Until bContinuelnner = True 'end of continue loop 
Next i 

dUBest = -dlnfinity 

For i = 1 To (iNumMOS - 1) 

If (dUArray(i) > dUBest) Then 
dUBest = dUArray(i) 
ilndex = i 
End If 
Next i 

bContinuelnner = True 
Do 'bContinuelnner do loop 
If (dUBest < (-dTest)) Then 
NewCostl = dOldCostl 
Exit Do 'exits bContinuelnner do loop 
End If 

NewCostl = dOldCostl - dCost(iIndex, dSoln(ilndex)) + dCost(iIndex, 
(dSoln(ilndex) + 1)) 
dSoln(ilndex) = dSoln(ilndex) + 1 
Loop Until bContinuelnner = True 
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End Function 


'WriteResults is a subroutine that will write the final solution to the Opt Solution 
'worksheet. 

'@param - rMOSData is a range variable that contains data for all the MOSZs. 

- rBudgetlnfo is a range variable that contains the budget and pay 
information needed for the SRB Optimization problem. 

- dSoln is an array of doubles that contains the multipliers solution set. 

Sub WriteResults(rMOSData As Range, rBudgetlnfo As Range, dSoln() As Double) 
Dim i As Integer 
Dim rOutput As Range 

Set rOutput = wsOptSoln.Range("OutputStart", wsOptSoln.Range("OutputStart") 
.Offset(iNumMOS, 10)) 
rOutput.Clear 
rOutput(l, 1) = "Occ Field" 
rOutput(l, 2) = "MOS" 
rOutput(l, 3) = "Zone" 
rOutput(l, 4) = "EASPop" 
rOutput(l, 5) = "B/Space" 
rOutput(l, 6) = "Multiplier" 
rOutput(l, 7) = "ExpReupNumber" 
rOutput(l, 8) = "RndReUpNumber" 
rOutput(l, 9) = "Over/Under" 
rOutput(l, 10) = "Expected Cost" 

For i = 2 To iNumMOS 

rOutput(i, 1) = rMOSData(i, 1) 
rOutput(i, 2) = rMOSData(i, 2) 
rOutput(i, 3) = rMOSData(i, 3) 
rOutput(i, 4) = rMOSData(i, 5) 
rOutput(i, 5) = rMOSData(i, 8) 
rOutput(i, 6) = (dSoln(i - 1) - 1) * dGap 

rOutput(i, 7) = rMOSData(i, (rOutput(i, 6) / dGap + 1) + 9) * rOutput(i, 4) 

Call FormatTwoDecimal(rOutput(i, 7)) 

rOutput(i, 8) = Application.WorksheetFunction.Round(rOutput(i, 7), 0) 
rOutput(i, 9) = rOutput(i, 8) - rOutput(i, 5) 

rOutput(i, 10) = rBudgetInfo(3, 1) * rBudgetInfo((rOutput(i, 3)) + 4, 1) * _ 
rOutput(i, 6) * rOutput(i, 7) 

Call FormatTwoDecimal(rOutput(i, 10)) 

Next i 
End Sub 

'File: Utilities Module in VBA 
'Date Created: 17 June 2007 
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'Last Updated: 25 July 2007 

t 

'Utilities is a Module that contains various formatting subroutines. 

f 

'@ Author Kent Robbins, Jr. 

t 

Option Explicit 

'FormatTwoDecimal is a subroutine that will format a given range to two decimal 
'places. 

t 

'@param - rTable is the range of values to be formatted. 

f 

Sub FormatTwoDecimal(rTable As Range) 
rTable. NumberFormat = "0.00" 

End Sub 

'FormatOccFielclText is a subroutine that will format a given range to text and 
'align it to the right of the cell. 

f 

'@param - rTable is the range of cells to be formatted. 

f 

Sub FormatOccFieldText(rTable As Range) 
rTable. NumberFormat = "@" 
rTable. HorizontalAlignment = xlRight 
End Sub 

File: WriteFiles Module in VBA 
'Date Created: 17 June 2007 
'Last Updated: 25 July 2007 

f 

'WriteFiles is a VBA Modula that will write will read in the GAMS solutions and 
'write the solutions along with information about the MOSs to the OptSoln worksheet. 

f 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 

WriteDataMaster is a subroutine that will gather information about MOSs that 
'is needed by the GAMS solver and write that informatino into .clat files. 

'@param - none 

f 

Sub WriteDataMaster() 

Dim i As Integer 
Dim j As Integer 
Dim rMOSData As Range 
Dim rAvgPays As Range 
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Dim rScalarData As Range 

Dim sFileNames() As String 

Dim slndex() As String 

Dim iLength As Integer 

Dim dMaxTrnCost As Double 

ChDrive Left$(ActiveWorkbook.Path, 1) 

ChDir ActiveWorkbook.Path 

Set rMOSData = wsReducedProblem.Range(wsReducedProblem.Range _ 
("RFPOutputStart").End(xlToRight), wsReducedProblem.Range _ 
("RFPOutputStart").End(xlDown)) 

Set rAvgPays = wsReducedProblem.Range("AvgPayA", wsReducedProblem _ 
.Range("AvgPayA").Offset(3, 0)) 

Set rScalarData = wsReducedProblem.Range("OverageWeight", wsReducedProblem 
.Range("OverageWeight").Offset(3, 0)) 
iLength = rMOSData.Rows.Count 
dMaxTrnCost = 0# 

ReDim sFileNames(l To 15) As String 
ReDim slndex(l To iLength) As String 
sFileNames(l) = "MOSZoneList.dat" 
sFileNames(2) = "MultiplierList.dat" 
sFileNames(3) = "AvgPay.dat" 
sFileNames(4) = "MOSPop.dat" 
sFileNames(5) = "EASPop.dat" 
sFileNames(6) = "PresetMult.dat" 
sFileNames(7) = "WgtFactor.dat" 
sFileNames(8) = "BoatSpace.dat" 
sFileNames(9) = "TrnCost.dat" 
sFileNames(lO) = "ReenlistmentRates.dat" 
sFileNames(l 1) = "MultiplierNumbers.dat" 
sFileNames(12) = "Q.dat" 
sFileNames(13) = "Budget.dat" 
sFileNames(14) = "Yrs.dat" 
sFileNames(15) = "MaxTrnCost.dat" 

For i = 2 To iLength 

If (rMOSData(i, 3) = 1) Then 

slndex(i) =.+ CStr(rMOSData(i, 2).Value) + "a" +. 

Elself (rMOSData(i, 3) = 2) Then 

slndex(i) =.+ CStr(rMOSData(i, 2).Value) + "b" +. 

Else 

slndex(i) = + CStr(rMOSData(i, 2).Value) + "c" + 

End If 

If (rMOSData(i, 9) > dMaxTrnCost) Then 
dMaxTrnCost = rMOSData(i, 9) 

End If 
Next i 
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For i = 1 To 11 

Call WriteData(i, sFileNames(i), slndex(), rMOSData, rAvgPays, iFength, 
rScalarData) 

Next i 

For i = 12 To 15 

Dim dSend As Double 
If i = 15 Then 

dSend = dMaxTrnCost 
Else 

dSend = rScalarData(i - 11) 

End If 

Call WriteDataScalar(sFileNames(i), dSend) 

Next i 
End Sub 

'WriteDataScalar is a subroutine that will write scalar data into a .dat file for 
'GAMS. 

t 

'@param - none 

f 

Sub WriteDataScalar(sNames As String, dData As Double) 

Dim iFN As Integer 
iFN = FreeFile() 

Open sNames For Output As iFN 
Print #iFN, "/" & dData & 7" 

Close iFN 
End Sub 

'WriteData is a subroutine that will write MOSZ data into a .dat file for 
GAMS. 

t 

'@param - iPosition As Integer is the index for the files to be written. 

- sNames As String is the file name for the data to be written. 

- slndex() As String is the name of the MOSZ to be written. 

- rMOSData As Range is the range of data to be printed. 

- rPay As Range variable of the avg pays to be printed. 

- iLength As Integer is the number of files to be printed. 

- rScalarData As Range variable of all the scalar values needed. 

t 

Sub WriteData(iPosition As Integer, sNames As String, slndex() As String, _ 
rMOSData As Range, rPay As Range, iFength As Integer, _ 
rScalarData As Range) 

Dim i As Integer, j As Integer, iFN As Integer, iNumMult As Integer 
iFN = FreeFile() 

Open sNames For Output As iFN 
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If (iPosition =1) Then 
Print #iFN, 7" 

For i = 2 To iLength 
Print #iFN, slndex(i); 

If i <> iLength Then 
Print #iFN, 

End If 
Next i 

Print #iFN, 7" 

Elself (iPosition = 2) Then 
Print #iFN, 7" 

For i = 1 To (rMOSData.Columns.Count - (9)) 
Print #iFN, i; 

If i <> (rMOSData.Columns.Count - (9)) Then 
Print #iFN, 

End If 
Next i 

Print #iFN, "/" 

Elself (iPosition = 3) Then 
Print #iFN, 7" 

For i = 2 To iLength 
Print #iFN, slndex(i) & 

If (rMOSData(i, 3) = 1) Then 
Print #iFN, rPay(l) 

Elself (rMOSData(i, 3) = 2) Then 
Print #iFN, rPay(2) 

Else 

Print #iFN, rPay(3) 

End If 

If (i <> iLength) Then 
Print #iFN, 

End If 
Next i 

Print #iFN, 7" 

Elself (iPosition =10) Then 

iNumMult = rMOSData.Columns.Count - 9 
Print #iFN, "dummy"; 

For i = 1 To iNumMult 
Print #iFN, & i; 

Next i 
Print #iFN, 

For i = 2 To iLength 
Print #iFN, slndex(i); 

For j = iPosition To rMOSData.Columns.Count 
Print #iFN, & rMOSData(i, j); 
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Next j 
Print #iFN, 

Next i 

Elself (iPosition =11) Then 

iNumMult = rMOSData.Columns.Count - 9 
Print #iFN, "/" 

For i = 1 To iNumMult 

Print #iFN, i & "=" & rMOSData(l, 9 + i); 

If (i <> iNumMult) Then 
Print #iFN, 

End If 
Next i 

Print #iFN, "/" 

Elself (iPosition =12) Then 
Print #iFN, "/" 

For i = 1 To 4 

Print #iFN, i & "=" & rScalarData(i); 

If (i <> 4) Then 
Print #iFN, 

End If 
Next i 

Print #iFN, "/" 

Else 

Print #iFN, "/" 

For i = 2 To iLength 

Print #iFN, slndex(i) & "=" & rMOSData(i, iPosition); 

If (i <> iLength) Then 
Print #iFN, 

End If 
Next i 

Print #iFN, "/" 

End If 
Close iFN 
End Sub 

'File: Solve Module in VBA 
'Date Created: 17 June 2007 
'Last Updated: 30 July 2007 

t 

'SolveGAMS is a subroutine that will run the GAMS solution model in a shell 
'program. 

r 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 
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Sub SolveGAMS() 

ChDrive Left$(ActiveWorkbook.Path, 1) 

ChDir ActiveWorkbook.Path 

Shell "gams ThesisSRBModel.gms", vbNormalFocus 
End Sub 

'File: ReadResults Module in VBA 
'Date Created: 17 June 2007 
'Last Updated: 14 July 2007 

f 

'ReadResults is a VBA Module that will read in the optimum SRB multiplier results from 
'GAMS and print them to Opt Solution worksheet along with its associated data. 

f 

'@ Author Kent Robbins, Jr. 

f 

Option Explicit 

'GAMSResults is a subroutine that will read in the optimum SRB multiplier results as 
'an outputfrom GAMS. Results will be outputted with other important data for the 
'MOS's. 

t 

'@param - None 

f 

Sub GAMSResults() 

Dim iFN As Integer 

Dim i As Integer 

Dim j As Integer 

Dim iLength As Integer 

Dim r As Double 

Dim dGap As Double 

Dim dExpReUpNum As Double 

Dim rMOSData As Range 

Dim rOutput As Range 

Dim rAvgPays As Range 

ChDrive Left$(ActiveWorkbook.Path, 1) 

ChDir ActiveWorkbook.Path 

Set rMOSData = wsReducedProblem.Range(wsReducedProblem.Range _ 
("RFPOutputStart").End(xlToRight), wsReducedProblem.Range _ 
("RFPOutputStart").Fnd(xlDown)) 

Set rAvgPays = wsReducedProblem.Range("AvgPayA", wsReducedProblem _ 
.Range("AvgPayA").Offset(3, 0)) 
iLength = rMOSData.Rows.Count 
dGap = rMOSData(l, 11) - rMOSData(l, 10) 

Set rOutput = wsOptSoln.Range("OutputStart", wsOptSoln.Range("OutputStart") _ 
.Offset(iLength, 10)) 
rOutput(l, 1) = "Occ Field" 
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rOutput(l, 2) = "MOS" 
rOutput(l, 3) = "Zone" 
rOutput(l, 4) = "EASPop" 
rOutput(l, 5) = "B/Space" 
rOutput(l, 6) = "Multiplier" 
rOutput(l, 7) = "ExpReupNumber" 
rOutput(l, 8) = "RndReUpNumber" 
rOutput(l, 9) = "Short/Over" 
rOutput(l, 10) = "Exact Expected Cost" 
iFN = FreeFile() 

Open "SRBOptResults.csv" For Input As iFN 
For i = 2 To iFength 

rOutput(i, 1) = rMOSData(i, 1) 
rOutput(i, 2) = rMOSData(i, 2) 
rOutput(i, 3) = rMOSData(i, 3) 
rOutput(i, 4) = rMOSData(i, 5) 
rOutput(i, 5) = rMOSData(i, 8) 

Input #iFN, r 

dExpReUpNum = rMOSData(i, (r / dGap + 1) + 9) * rOutput(i, 4) 

If dExpReUpNum = 0 Then 
rOutput(i, 6) = 0 

dExpReUpNum = rMOSData(i, 10) * rOutput(i, 4) 

Else 

rOutput(i, 6) = r 
End If 

rOutput(i, 7) = dExpReUpNum 
Call FormatTwoDecimal(rOutput(i, 7)) 

rOutput(i, 8) = Application.WorksheetFunction.Round(rOutput(i, 7), 0) 
rOutput(i, 9) = rOutput(i, 8) - rOutput(i, 5) 

rOutput(i, 10) = 4 * rAvgPays(rOutput(i, 3)) * rOutput(i, 6) * rOutput(i, 7) 
Call FormatTwoDecimal(rOutput(i, 10)) 

Input #iFN, r 
Next i 
Close iFN 
End Sub 
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APPENDIX C: EXAMPLE INPUTS AND OUTPUTS 


This appendix contains the examples of the inputs and outputs of the Excel 
interface that was used to solve the SRB multiplier selection problem for FY04. 
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Figure 1. Excel Worksheet for Zone A Reenlistment Data. 
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Figure 2. Excel Worksheet of complied MOSZs and data 
need for the execution of the Model. 
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Figure 3. Excel Worksheet with MOSZ data, multiplier 

solution and expected cost. 
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